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ABSTRACT 


Possible  improvements  to  the  minimum  resolvable  temperature  difference  (MRTD) 
entered  into  TAWS  are  considered.  FLIR92  is  modified  to  include  atmospheric  turbu¬ 
lence  that  depends  on  height  in  the  atmosphere.  Resultant  MRTDs  are  compared  to 
the  operational  FLIR92  MRTD  predictions  excluding  atmospheric  turbulence.  The 
difference  in  the  MRTD  is  only  apparent  in  the  higher  frequency  regime  and  is  less 
than  0.05%  of  the  MRTD  values  for  a  typical  test  case.  MRTD  is  calculated  by 
FLIR92  and  NVThermlP  over  desert  and  marine  locations  and  the  resultant  MRTDs 
entered  into  TAWS  to  compare  maximum  detection  range.  NVThermlP  yielded  a 
larger  maximum  detection  range  by  up  to  1.5%  over  the  desert  and  2%  over  water. 
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I. 


INTRODUCTION 


Target  Acquisition  Weapons  Software  (TAWS)  is  used  operationally  to  predict 
the  performance  of  electro-optical  weapons  and  their  navigation  systems.  For  sensors 
operating  in  the  infrared  (IR)  spectral  band,  the  minimum  resolvable  temperature  dif¬ 
ference  (MRTD)  is  determined  from  a  desktop  computer  model  such  as  the  Forward 
Looking  Infrared  92  (FLIR92)  or  the  Night  Vision  Thermal  and  Image  Processing 
model  (NVThermlP).  Both  models  are  capable  of  providing  an  MRTD,  which  may 
be  manipulated  for  use  in  TAWS,  but  there  are  differences  in  their  calculations  that 
become  apparent  at  high  and  low  frequencies.  Also,  although  FLIR92  does  not  oper¬ 
ationally  include  the  effects  of  atmospheric  turbulence,  NVThermlP  does  for  limited 
conditions. 

In  this  study,  FLIR92  and  NVThermlP  are  compared  in  their  current  opera¬ 
tional  forms  for  various  atmospheric  conditions.  Then  FLIR92  is  modified  to  include 
atmospheric  turbulence  and  the  two  version  of  FLIR92  are  compared  for  the  varied 
atmospheric  conditions. 

All  acronyms  in  this  report  are  listed  in  alphabetical  order  in  Appendix  E. 

A.  MOTIVATION  FOR  COMPARISON 
1.  Current  Status  of  TAWS 

Currently,  the  MRTD  used  in  TAWS  predictions  of  target  detection  by  IR 
sensor  is  calculated  by  FLIR92  from  basic  system  parameters.  While  this  has  been 
sufficient  for  previous  iterations  of  TAWS,  recent  improvements  in  model  resolution 
and  the  image  resolution  of  sensors  being  modeled  suggest  that  perhaps  NVThermlP 
may  be  a  better  choice.  NVThermlP  incorporates  a  number  of  improvements  intended 
to  improve  the  determined  MRTDs  for  the  sensors  in  the  extrema  of  low  and  high 
spatial  frequencies,  in  other  words  for  very  large  and  very  small  targets.  NVThermlP 
also  includes  optical  turbulence,  which  should  improve  the  MRTD  calculation.  With 
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increased  sensitivities  of  sensors,  including  optical  turbulence  should  improve  the 
MRTD  prediction.  Since  FLIR92  and  NVThermlP  behave  differently  at  high  and  low 
frequencies  in  marine  or  desert  environments,  one  model  may  substantially  improve 
the  MRTD  calculation,  and  thus  the  TAWS  predictions,  over  the  other  model.  In 
particular,  the  temperature  difference  between  the  surface  and  the  ambient  air  may 
influence  model  calculations  and  this  temperature  difference  is  much  greater  during 
the  day  over  a  desert. 

2.  Differences  Between  FLIR92  and  NVThermlP 

Both  FLIR92  and  NVThermlP  are  desktop  computer  models  developed  to 
predict  standard  summary  performance  figures  of  merit  for  thermal  imaging  sys¬ 
tems.  They  are  PC  based  programs  that  model  passive  sensors  which  detect  emitted 
and  reflected  radiation.  Using  basic  system-level  parameters,  these  models  calculate 
the  modulation  transfer  function  (MTF),  the  noise  equivalent  temperature  difference 
(NETD),  the  MRTD,  and  the  minimum  detectable  temperature  difference  (MDTD) 
for  sensors  looking  at  specific  targets.  FLIR92  and  NVThermlP  both  predict  the 
MRTD  that  a  human  can  discriminate  when  using  a  thermal  imager  operationally. 
NVThermlP  also  predicts  the  range  at  which  target  acquisition  will  be  successful, 
using  the  specified  thermal  imager.  Thus,  the  output  from  each  model  is  used  to 
determine  whether  or  not  a  system  will  achieve  the  MTF,  system  noise,  MRTD,  and 
MDTD  that  is  required  to  effectively  perform  a  given  mission.  For  each  mission,  the 
conditions  necessary  to  meet  the  target  acquisition  and  discrimination  requirements 
may  vary. 

FLIR92  and  NVThermlP  share  certain  basic  assumptions.  These  assumptions 
are  stated  in  the  NVThermlP  User’s  Manual  (NVESD  2001),  and  follow  the  steps  of 
modeling  FLIR  systems  outlined  in  the  Infrared  and  Electro- Optical  (IREO)  Hand¬ 
book  Vol.4  (1993).  All  MTFs  are  assumed  to  be  separable,  so  the  total  system  MTF 
is  calculated  as  the  product  of  all  sub-system  MTFs.  This  approach  reduces  calcula¬ 
tion  complexities  because  the  analysis  is  simplified  to  one  dimension  and  cross-terms 
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are  eliminated.  However,  there  is  almost  always  some  calculation  error  associated 
with  assuming  separability.  This  error  becomes  significant  in  certain  cases  such  as  for 
diamond-shaped  detectors,  but  is  neglected  in  all  FLIR92  and  NVThermlP  calcula¬ 
tions  regardless  of  detector  shape  and  type  (NVESD  2001). 

An  MTF  describes  the  spatial  frequency  response  of  a  system.  It  is  the  contrast 
at  a  given  spatial  frequency  compared  to  the  contrast  at  a  lower  frequency.  The  total 
system  MTF  is  a  product  of  all  component  MTFs.  Each  MTF  may  be  calculated 
as  the  Fourier  transform  of  the  point  spread  function  (PSF)  of  the  component,  or 
the  response  of  the  imaging  component  to  a  point  of  light.  An  imaging  system 
will  experience  some  degradation  of  the  image  due  to  imperfections  in  the  optics, 
electronics,  or  even  the  observer  eye  and  this  degradation  is  what  the  MTF  describes. 
If  some  system  component  did  not  contribute  to  the  degradation,  then  the  MTF  for 
that  component  would  be  unity  (Driggers  1999). 

Both  FLIR92  and  NVThermlP  assume  all  blurs  to  be  symmetrical  in  order  to 
keep  all  MTFs  real.  The  models  assume  that  there  exists  a  region  of  the  field  of  view 
(FOV)  that  is  isoplanatic.  Goodman  (1968)  explains  that  in  an  isoplanatic,  or  space- 
invariant,  linear  imaging  system  the  image  of  a  point-source  will  change  location, 
but  not  functional  form,  as  the  point  source  moves.  The  approximately  isoplanatic 
region  of  the  FOV  is  modeled  in  FLIR92  and  NVThermlP  by  a  linear  shift-invariant 
process  and  the  MTF  in  that  region  is  approximated  by  symmetrical  blur.  The  blur 
is  the  result  of  real  world  effects  on  the  image  and  can  be  clue  to  aberrations  and 
other  manufacturing  defects  in  the  optical  system.  For  a  point  source,  the  image 
is  called  a  blur  circle  because  diffraction,  aberrations,  manufacturing  defects,  and 
assembly  tolerances  prevent  perfect  resolution  of  the  singularity  (Driggers  1999).  In 
FLIR92,  the  blur  MTF  is  geometric  and  is  discussed  in  further  detail  in  the  next 
chapter.  The  symmetrical  blur  approximation  is  not  an  accurate  reflection  of  real 
world  systems,  so  the  blur  approximation  for  the  optics  is  not  completely  correct. 
Similarly,  assuming  symmetrical  blurs  does  not  accurately  reflect  the  electronics  used 
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in  the  various  sensor  systems  under  consideration:  a  low  pass  filter  would  not  result 
in  a  time  delay  or  phase  shift  using  this  assumption,  for  example  (NVESD  2001). 

NVThcrmlP  provides  some  new  capabilities  beyond  FLIR92.  NVThermlP 
provides  target  acquisition  performance  predictions  for  staring  imagers,  not  just  the 
first  and  second  generation  thermal  scanning  sensors.  In  the  NVThcrmlP  calculation, 
the  MTF  representing  the  function  of  the  human  eye  includes  factors  that  are  ignored 
in  FLIR92.  The  MRTD  prediction  is  changed  from  FLIR92;  NVThermlP  produces 
contrast  transfer  functions  (CTFs)  as  the  primary  output.  A  laboratory  MRTD  is  cal¬ 
culated,  but  is  not  comparable  to  the  MRTD  output  from  FLIR92.  Also,  NVThermlP 
now  uses  the  Targeting  Task  Performance  (TTP)  metric  to  predict  the  probability 
of  target  detection,  recognition,  and  identification  at  given  ranges.  FLIR92  uses  the 
Johnson  criteria,  but  does  not  predict  probabilities  of  detection,  recognition,  or  iden¬ 
tification  (NVESD  2001).  In  this  report,  the  TTP  capabilities  in  NVThermlP  are 
neglected  since  only  the  MRTD  results  from  FLIR92  are  compared  to  those  from 
NVThermlP. 

3.  Advantages  and  Disadvantages  of  FLIR92 

a.  Advantages  of  FLIR92 

Since  NVThermlP  uses  the  same  basic  MRTD  prediction  theory  with 
modifications  to  improve  on  FLIR92  (NVESD  2001),  familiarity  would  be  the  primary 
benefit  of  continuing  with  FLIR92.  These  modifications  are  discussed  in  the  next 
sections,  regarding  the  advantages  of  NVThermlP.  If  FLIR92  operates  at  a  sufficiently 
accurate  and  precise  level  to  satisfy  the  increasingly  more  stringent  requirements  in 
operational  use  of  TAWS,  then  there  is  no  reason  to  change.  Also,  NVThermlP 
does  not  calculate  a  MDTD,  which  is  more  useful  for  indicating  thermal  imager 
performance  for  point  sources  and  aperiodic  targets.  TAWS  uses  the  MDTD  for 
certain  calculations,  when  available.  For  these  two  reasons,  if  a  simple  modification 
to  FLIR92  to  include  optical  turbulence  improves  the  FLIR92  MRTD  prediction,  then 
it  would  be  advantageous  to  continue  using  FLIR92. 
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b.  Disadvantages  of  FLIR92 

There  are  two  major  characteristics  of  FLIR92  that  may  lead  to  errors 
in  performance  predictions  for  staring  arrays  in  sensors  being  modeled.  Improved 
sensitivities  in  staring  arrays  mean  that  the  limitations  of  the  eye  in  detecting  con¬ 
trasts  is  now  restricting  the  performance  of  these  sensors  (NVESD  2001),  but  FL1R92 
MRTD  calculations  do  not  include  any  adjustments  due  to  less-than-perfect  contrast 
perceptions  of  the  eye.  Likewise,  limitations  on  detector  size,  spacing,  and  fill  factor 
(ratio  of  active  cell  area  to  total  array  area  (Driggers  1999))  can  cause  under-sampled 
imagery  for  staring  sensors.  To  avoid  the  under-sampled  imagery,  FL1R92  has  an 
absolute  cutoff  at  the  half  sample  rates  of  imagers  (Nyquist  frequency),  but  when 
used  with  the  Johnson  criteria,  it  can  lead  to  incorrectly  negative  predictions  for 
most  staring  imagers  (NVESD  2001).  FL1R92  also  does  not  account  for  atmospheric 
turbulence,  which  may  be  particularly  significant  in  the  highly  turbulent  regime  just 
above  a  hot  noon-time  desert. 

4.  Advantages  and  Disadvantages  of  NVThermlP 

a.  Advantages  of  NVThermlP 

NVThermlP  is  considered  an  advance  over  FL1R92  primarily  because  it 
addresses  the  two  possible  errors  in  FL1R92  performance  predictions  regarding  staring 
arrays,  as  discussed  above.  NVThermlP  also  does  include  atmospheric  turbulence  in 
a  fairly  rudimentary  form. 

b.  Disadvantages  of  NVThermlP 

There  are  three  major  hesitations  in  switching  to  use  NVThermlP  with 
TAWS.  First  is  a  lack  of  familiarity.  Otherwise,  since  NVThermlP  was  developed  to 
address  weaknesses  in  FL1R92,  it  should  be  advantageous  to  use.  For  instance,  NV¬ 
ThermlP  uses  the  Targeted  Task  Performance  (TTP)  Metric  instead  of  the  Johnson 
Metric  to  provide  a  better  performance  estimate.  This  result  is  then  used,  like  with 
FL1R92,  to  predict  the  target  acquisition  range  performance  for  given  sensors.  Sec¬ 
ond,  while  NVThermlP  predicts  the  MRTD,  it  is  a  laboratory  MRTD  that  is  not 
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comparable  to  the  MRTD  from  FLIR92.  Instead,  NVThermlP  predicts  the  CTF 
(NVESD  2001),  which  must  be  manipulated  before  insertion  into  TAWS.  Finally, 
NVTHermlP  does  not  produce  an  MDTD,  but  TAWS  uses  the  MDTD  for  certain 
applications. 

5.  Phenomena  Missing  from  FLIR92  and  NVThermlP 

Despite  improvements  in  NVThermlP,  there  are  details  missing  from  both 
FLIR92  and  NVThermlP.  One  important  lack  is  the  optical  turbulence.  NVThermlP 
does  include  an  average  optical  transmission  input,  but  neither  model  considers  the 
effect  of  optical  turbulence  variations  along  the  optical  path  from  the  target  to  the 
sensor.  The  next  section  discusses  this  in  more  detail. 

Neither  FLIR92  nor  NVThermlP  have  any  adjustments  for  polarization  of  the 
target.  Since  the  system  response  varies  depending  on  the  polarization,  this  lack  will 
be  especially  apparent  as  imager  sensitivities  increase.  Both  models  work  best  for 
nearly  symmetric  targets,  but  may  poorly  represent  more  realistic  cases. 

B.  MOTIVATION  TO  TEST  IMPORTANCE  OF  OPTI¬ 
CAL  TURBULENCE 

Recently,  sensors  with  higher  resolution  have  been  developed.  These  IR  im¬ 
agers  are  more  sensitive  to  smaller  targets,  or  a  higher  spatial  frequency.  In  order 
to  effectively  model  the  responses  of  these  sensors,  it  is  becoming  necessary  to  con¬ 
sider  and  perhaps  include  previously  negligible  effects  like  optical  turbulence.  At  the 
sensitivities  of  the  new  sensors,  the  range  for  minimum  detectability  may  be  signif¬ 
icantly  affected  by  slight  variations  in  the  optical  turbulence  along  the  path  from 
the  target  to  the  sensor,  perhaps  even  as  regards  a  vertical  pathlength  through  the 
atmosphere.  Optical  turbulence  causes  blurring  around  the  edges  of  images  (Fante 
1980)  and  makes  it  more  difficult  to  resolve  the  image  of  the  target  well  enough  to 
identify  it.  According  to  the  IREO  Handbook  Vol.3  (1993),  the  effects  of  pressure 
variations  due  to  atmospheric  turbulence  are  typically  negligible  compared  to  those 
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due  to  temperature  and  humidity  fluctuations.  In  either  very  hot  atmospheres  such  as 
those  that  may  be  operationally  important  in  South  Asia,  or  in  very  sensitive  appli¬ 
cations  -  such  as  detection  of  a  very  small  target,  the  effect  of  atmospheric  turbulence 
may  become  significant.  Recent  improvements  in  sensor  responsivities  have  increased 
the  number  of  operations  in  the  sensitive  applications  category  and  therefore  exam¬ 
ining  the  potential  impact  of  optical  turbulence  along  the  optical  path  through  the 
atmosphere  is  crucial. 

C.  PROJECT  OVERVIEW 

Two  comparisons  to  test  possible  improvements  to  the  MRTD  calculation  were 
accomplished. 

First,  a  MatLab  version  of  the  operationally  used  FLIR92  was  modified  to  in¬ 
clude  variations  in  optical  turbulence  along  the  pathlength.  Then  the  original  FLIR92 
was  compared  to  the  modified  FLIR92  to  investigate  the  significance  of  improvements 
in  the  MRTD.  Resultant  MRTDs  were  inserted  into  TAWS  and  the  maximum  detec¬ 
tion  ranges  were  compared.  The  generic  IR  sensor  in  TAWS  was  used  and  marine 
and  desert  environments  were  tested. 

Second,  the  MRTDs  calculated  by  FLIR92  and  by  NVThermlP  were  com¬ 
pared.  All  input  parameters  were  held  constant  for  each  model.  The  resultant  MRTDs 
from  FLIR92  and  NVThermlP  were  entered  into  TAWS.  The  atmospheric  conditions 
were  varied  to  test  the  marine  and  desert  environments.  In  TAWS,  the  generic  IR 
sensor  was  used  for  all  tests,  modified  to  include  the  MRTD  from  either  FLIR92  or 
NVThermlP. 

1.  Comparison  Parameters 

To  compare  FLIR92  and  NVThermlP  under  the  conditions  described,  param¬ 
eters  calculated  by  each  were  compared.  Factors  to  compare  and  evaluate  in  order  to 
determine  the  significance  of  any  differences  are  the  horizontal  and  vertical  MRTDs 
and  system  MTFs  at  specified  spatial  frequencies.  Again,  the  extrema  of  frequencies 
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are  the  regions  of  particular  interest  in  this  study.  Because  NVThermlP  outputs  the 
CTF,  in  order  to  compare,  the  resultant  CTF  was  multiplied  by  2  times  the  scene 
contrast  temperature  (SCT)  to  obtain  the  MRTD  (Jon  Hixson,  personal  communica¬ 
tion). 

The  terms  ‘horizontal’  and  ‘vertical’  directions  as  used  in  this  report  refer  to 
the  along-bar  and  cross-bar  directions  of  a  four-bar  target,  as  shown  in  Figure  1. 


h  axis 

Horizontal  Target 
Orientation 


Vertical  Target 
Orientation 


Figure  1.  Horizontal  and  Vertical  Directions  with  Respect  to  a  4-Bar  Target. 


2.  Comparisons  in  TAWS 

Once  the  MRTDs  were  determined  from  FLIR92  and  NVThermlP,  the  re¬ 
sults  were  entered  into  TAWS.  Then,  the  range  output  from  TAWS  was  compared 
for  the  various  different  iterations  to  determine  how  much  and  how  significant  the 
improvements  to  the  maximum  detection  range  in  TAWS  were. 
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II.  THERMAL  IMAGING  SYSTEMS 
PERFORMANCE  MODELS 

A.  FLIR92 

1.  MTFs  Included  in  FLIR92 

FLIR92  uses  a  total  system  MTF  to  calculate  the  MRTD.  The  total  system 
MTF  is  the  product  of  the  prefilter  and  postfilter  MTFs,  including  the  noise  filter 
MTFs.  For  most  component  MTFs,  FLIR92  can  supply  a  model  prediction,  but 
allows  for  a  user  to  directly  input  a  measured  or  previously  known  MTF  for  the 
component.  Generally,  the  measured  MTFs  improve  performance  estimates  and  are 
preferred  over  the  MTFs  predicted  and  calculated  by  FLIR92. 

In  FLIR92,  the  horizontal  and  vertical  prefilter  MTFs  for  a  sensor  are  the 
products  of  the  optical,  detector  spatial,  focal  plane  array  (FPA)  integration,  and 
image  motion  MTFs  in  the  horizontal  and  vertical  directions  respectively.  FLIR92 
allows  user-defined  MTFs  for  all  but  the  FPA  integration  MTF. 

For  scanned  observing  systems,  filtering  in  the  temporal  domain  must  be  trans¬ 
lated  into  the  spatial  domain.  The  detector  temporal  MTF,  electronics  high  and  low 
frequency  response  MTFs,  and  boosting  MTF  are  included  in  the  temporal  postfilter 
MTF  that  describes  this  translation.  The  temporal  postfilter  MTF  in  FLIR92  only 
applies  for  scanning  systems  and  is  only  defined  in  the  horizontal  direction.  FLIR92 
allows  a  user-defined  temporal  postfilter  MTF. 

Both  scanning  and  staring  systems  have  spatial  postfilters  that  model  the 
system  processing  components.  The  spatial  postfilter  MTFs  are  the  electro-optical 
multiplexer  MTF,  the  digital  filter  MTF,  the  display  MTF,  the  Charge-Coupled  De¬ 
vice  (CCD)  charge  transfer  efficiency  MTF,  the  display  sample  and  hold  MTF,  and 
the  eye  MTF.  All  are  defined  in  both  the  vertical  and  horizontal  direction. 

FLIR92  allows  for  an  extra  user-determined  prefilter  MTF  and  three  extra 
postfilter  MTFs  to  enable  adequate  description  of  the  system.  All  the  following 
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MTFs  are  standard  in  FLIR92  and  are  summarized  from  the  Analyst’s  Reference 
Guide  (1993)  and  compared  against  original  sources  or  the  IREO  Handbook  (1993) 
as  specified.  For  scanning  systems,  the  scan  direction  is  assumed  to  be  horizontal  to 
the  bars  of  the  target. 

a.  Optical  MTF 

The  optical  MTF  includes  the  diffraction-limited  MTF  and  geometric 
blur  MTF.  When  available,  the  optical  MTF  may  be  replaced  by  a  user-defined 
version  determined  from  direct  measurements  or  ray  tracing  program  predictions. 
The  diffraction-limited  MTF  describes  the  resolution  limitations  due  to  diffraction  in 
the  optics  of  the  sensor  and  is  valid  for  a  system  with  a  circular,  clear  aperture  ( IREO 
Handbook  Vol.4  1993).  It  it  of  the  form 


HMs) 


2 

71 


/ 

arccos 

V 


(  A/s 
VA) 


(2.1) 


where  A  is  the  wavelength  for  diffraction  in  /irn  and  D o  is  the  optics  aperture  diameter 
in  mm.  Since  arccos(a;)  is  only  physical  when  —  1  <  x  <  1,  then  the  spatial  frequency 
fs,  in  cycles/mrad,  must  satisfy  fs  <  A.  The  diffraction  limited  MTF  in  the  hori¬ 
zontal  direction  is  illustrated  in  Figure  2.  For  this  report,  the  same  inputs  are  used 
in  the  horizontal  and  vertical  directions,  so  the  MTF  in  the  vertical  direction  is  the 
same. 

Apertures  that  are  partially  obscured  or  are  not  perfectly  circular,  or 
aberrations  in  the  optics  may  cause  blurring  of  the  image  (Goodman  1968).  The 
optics  MTF  includes  a  geometric  blur  term  to  describe  this  phenomenon.  If  there  is 
no  user  input  for  the  geometric  blur  MTF,  then  FLIR92  assumes  the  blur  is  Gaussian 
and  the  MTF  has  the  form 


Hgb{fs)  =  exp(-27rV/2). 


(2.2a) 
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a,  in  mrad,  is  the  standard  deviation  for  a  circular  Gaussian  blur  distribution.  The 
IREO  Handbook  Vol.4  (1993)  clarifies:  a2  is 


cr 


2 


(2.2b) 


where  w  is  the  Gaussian  blur  spot  diameter  in  mrad  at  the  \  point. 

Figure  3  shows  the  geometric  blur  MTF  in  the  horizontal  direction. 
Again,  with  the  inputs  used  in  this  report,  the  MTF  is  identical  in  the  vertical  direc¬ 
tion. 

b.  Detector  Spatial  MTF 

Also  included  in  the  prefilter  MTF  is  the  detector  spatial  MTF  that 
compensates  for  the  finite  size  of  the  detector.  FL1R92  assumes  a  rectangular  detec¬ 
tor  geometry.  Like  for  the  optical  MTF,  the  IREO  Handbook  Vol.5  emphasizes  the 
importance  of  a  user-specified  detector  spatial  MTF,  particularly  if  the  rectangular 
detector  geometry  approximation  is  grossly  inaccurate. 


Diffraction  MTF  H 


Figure  2.  Diffraction  Limited  MTF  in  the  Horizontal  Direction. 
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A  Signal  Processing  in  the  Element  (SPRITE)  system,  incorporating 
the  signal  time  delay  and  integration  into  one  elongated  detector  element  ( IREO 
Handbook  Vol.3  1993)  has  a  horizontal  detector  spatial  MTF  of  the  form 

sm(nlhfs)  (  1  \ 


Hds(fs )  = 


irlhfs  ll  +  (2v rQfs)2 


(2.3a) 


where  Q  is  the  diffusion  length  and  A  is  the  instantaneous  field  of  view  (IFOV)  of 
the  readout  length,  both  in  mrad.  The  IFOV  is  the  solid  angle  through  which  a 
detector  is  sensitive  to  radiation  and  in  a  scanning  system.  This  is  the  subtended 
solid  angle  when  all  scanning  motion  stops.  The  vertical  detector  spatial  MTF  for 
SPRITE  systems  and  both  horizontal  and  vertical  MTFs  for  non-SPRITE  scanning 
and  staring  systems  is 


Hds(fs)  = 


sin(7 T6zfs) 

7T  Szfs 


(2.3b) 


where  Sz  is  the  detector  IFOV  in  mrad  (Boreman  1990).  The  IREO  Handbook  Vol.4 
(1993)  also  points  out  that  this  form  is  only  valid  when  the  point  spread  function  of 
the  detector  is  a  rectangular  function  in  direct  space,  but  it  is  a  reasonable  approxi- 


Blur  MTF  H 


Figure  3.  Geometric  Blur  MTF  in  the  Horizontal  Direction. 
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mation  in  most  cases.  The  MatLab  version  of  FLIR92  used  in  this  analysis  assumes 
a  non-SPRITE  system  and  uses  Equation  2.3b.  Figure  4  shows  the  form  of  the  de¬ 
tector  spatial  MTF  in  the  horizontal  direction,  and  again  it  is  identical  in  the  vertical 
direction. 


Det  aprtr  MTF  H 


Figure  4.  Detector  Spatial  MTF  in  the  Horizontal  Direction. 
This  plot  is  assuming  a  non-SPRITE  scanning  or  staring  system. 


c.  Focal  Plane  Array  Integration  MTF 

For  scanning  systems,  the  prefilter  MTF  also  includes  a  focal  plane 
array  integration  MTF  introduced  in  the  horizontal  direction  by  the  finite  integration 
time  of  the  detector.  If  tt  is  the  detector  integration  time  in  s  and  vs  is  the  scan 
velocity  in  mrad/s,  the  MTF  has  the  form  of  diffraction,  or 


Hdi{fs) 


sin(7r/sustj) 

-/,  ^  s  ^  / 


(2.4) 


This  MTF  is  only  significant  for  scanning  systems  because  staring  systems  do  not 
require  compensation  for  finite  integration  times  since  they  do  not  move  during  de¬ 
tection. 
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d.  Sample-scene  Phase  MTF 

When  considering  sampled  systems,  the  target  location  on  the  sampling 
grid  can  cause  a  phase  effect,  or  aliasing  effect.  The  aliasing  of  high  frequency  noise 
into  the  bandpass  of  the  sampled  signal  can  produce  false  signals  or  greatly  increased 
noise  or  interference  ( IREO  Handbook  Vol.3  1993).  The  MTF  representing  this  phase 
effect  is 


Hssp{fs 


(2.5a) 


Here,  tn  is  the  Nyquist  frequency  in  cycles/mrad  and  ©z  is  the  phase  angle  (rad)  be¬ 
tween  the  target  under  consideration  and  the  detectors  operating  at  rjv-  The  Nyquist 
frequency,  or  Nyquist  limit,  is  the  frequency  at  which  no  useful  information  is  trans¬ 
mitted  and  is  taken  as  half  the  scene  sample  frequency.  It  is  given  as 


(2.5b) 


where  sz  is  the  number  of  samples  per  detector  IFOV.  Input  frequencies  above  the 
Nyquist  limit  are  likely  to  appear  as  aliasing  signals  at  the  lower  frequencies  ( IREO 
Handbook).  At  optimal  operation,  ©^  is  set  equal  to  0,  but  for  average  conditions  0Z 
is  set  equal  to  0.785  rad  (45°).  In  the  inputs  used  for  the  MatLab  version  of  FLIR92, 
©^  is  set  to  0,  or  for  optimal  operation.  This  results  in  an  MTF  in  the  horizontal  and 
the  vertical  direction  that  is  constant  at  unity  for  all  frequencies  since  cos(0)  is  unity. 

e.  Image  Motion  MTF 

Image  motion  MTFs  are  included  in  FLIR92  to  account  for  how  the 
thermal  system  moves  with  respect  to  the  scene  being  imaged.  Two  causes  of  these 
movements  may  be  linear  motion  of  the  sensor  system  or  vibration  of  the  sensor 
platform.  These  movements  are  represented  in  the  linear  image  motion  MTF,  random 
image  motion  MTF,  and  sinusoidal  image  motion  MTF. 

The  linear  motion  MTF  explains  smearing  of  the  scene  across  the  detec¬ 
tors  due  to  movement  of  the  system  platform.  How  significant  the  smear  is  depends 
on  how  fast  the  platform  is  moving  and  how  long  the  detectors  are  exposed  to  the 
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scene  during  platform  movement.  Given  vr  is  the  relative  image  velocity  in  mrad/s, 
the  linear  image  motion  MTF  may  be  described  by 

sm(nvr  tifs) 


Hml(fs)  = 


(2.6) 


7TVrtifs 

When  ti  is  not  specihed,  the  integration  time  defaults  to  the  held  time.  Typically,  the 
linear  image  motion  is  only  significant  for  staring  systems  because  scanning  systems 
expect  the  sensor  to  be  moving,  or  scanning,  during  detection.  The  MatLab  version 
of  FLIR92  sets  this  MTF  to  unity  for  all  frequencies. 

For  simplicity,  the  random  image  motion  is  assumed  to  be  Gaussian. 
The  IREO  Handbook  Vol.4  (1993)  notes  that  this  Gaussian  term  includes  the  random 
vibration  of  the  sensor  with  respect  to  the  scene,  uncertainty  in  the  sensor  position, 
and  other  causes  of  image  jitter.  The  associated  MTF  is 

Hmr(fs)  =  exp(-2  7r2aloSfs)  (2.7) 


where  cr2LOS  is  the  standard  deviation  of  the  total  random  vibration  from  all  sources, 
in  mrad.  Figure  5  illustrates  the  form  of  the  random  motion  MTF  in  the  horizontal 
direction.  The  vertical  direction  is  identical  with  the  given  inputs. 

Finally,  the  sinusoidal  vibration  of  the  platform  may  be  described  by  a 
zero-order  Bessel  function,  as  shown  in  the  IREO  Handbook  Vol.5  (1993), 


Hms(fs)  =  J0(2vr  Afs). 


(2.8) 


A  is  the  amplitude  of  the  assumed  sinusoidal  vibration  in  mrad.  The  zero-order  Bessel 
function  has  the  form  shown  in  Fig  6,  similar  to  a  damped  cosine  function.  Using  the 
given  inputs,  though,  A  is  set  to  zero  and  Hms  is  thus  held  constant  at  unity  for  all 

frequencies  since  Jo(0)  =  1. 

/.  Detector  Temporal  MTF 

Because  the  detector  has  a  finite  response  time,  the  system  has  a  low- 
pass  temporal  filter  given  by 

(2,9) 


HMt)  = 


1+(£, 
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where  ft  is  the  temporal  frequency  in  Hz  and  is  related  to  the  spatial  frequency  by 
the  scan  velocity  since  ft  =  fsvs.  fdt  is  the  3dB  frequency  in  Hz  of  the  detector 
temporal  response.  The  IREO  Handbook  Vol.4  further  elaborates  on  Equation  2.9  - 
the  finite  temporal  response  to  radiation  incident  on  the  detectors  causes  blurring  in 
the  scanning  direction,  which  is  why  H^t  is  included  in  the  horizontal  for  scanning 
systems. 

g.  Electronics  Frequency  Responses 

The  frequency  responses  for  systems  that  are  AC-coupled  are  given  by 
multiple  pole  RC  filters.  The  IREO  Handbook  (1993)  gives  the  single-RC  low-pass 
and  high-pass  filter  MTFs,  which  are  then  raised  to  the  nth  power,  where  n  is  the 
number  of  filter  poles,  to  arrive  at  Equations  2.10a  and  2.10b. 

For  the  low  frequency  response,  the  high  pass  filter  MTF  is 


Hehp(ft) 


(2.10a) 


Motion  Rnd  MTF  H 


Figure  5.  Random  Motion  MTF  in  the  Horizontal  Direction. 
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where  fehp  is  the  electronics  3dB  frequency  in  Hz.  The  low  frequency  response  is  only 
significant  at  extremely  low  frequencies,  so  FLIR92  does  not  include  it  in  the  total 
MTF,  but  Hehp  is  used  to  calculate  noise  bandwidths. 


For  the  high  frequency  response,  the  low  pass  filter  MTF  is 


(2.10b) 


where  like  in  the  high  pass  filter  case,  feip  is  the  3dB  frequency  in  Hz  and  n  is  the 
number  of  filter  poles.  This  high  frequency  response  MTF  is  included  in  the  total 
MTF  since  the  high  frequency  response  is  significant  at  a  greater  range  of  frequencies. 

h.  Boosting  MTF 

For  the  scanning  system,  an  aperture  correction  MTF,  or  electronic 
boosting  MTF,  is  required.  The  boost  emphasizes  higher  frequencies  to  compensate 
for  the  reduced  depth  of  modulation  that  typically  occurs  at  higher  frequencies  due  to 
the  less-than- ideal  aperture  response  ( IREO  Handbook  Vol.4  1993).  It  has  the  form 


(2.11) 
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Figure  6.  Zeroth  Order  Bessel  Function 
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fb  is  the  frequency  in  Hz  of  the  maximum  boost  and  Ba  is  the  amplitude  of  the  boost 
at  a  maximum  frequency  fmax ■  As  the  IREO  Handbook  Vol.4  observes,  the  boost 
MTF  is  applied  along  the  scan  direction. 

i.  Electro- Optical  Multiplexer  MTF 

For  the  LED  electro-optical  multiplexer,  the  MTF  is  of  the  form  of 


diffraction, 


Heom(fs) 


sin  {'Kfs8LED) 


(2.1.2) 

e0m  is  not  included 


f sd  LED 

and  Sled  is  the  angular  subtense  of  the  LED  element  in  mrad.  H, 
in  the  MatLab  version  of  FLIR92. 

j.  Digital  Filter  MTF 

Digital  Liters  have  a  linear  phase  symmetrical  impulse  response  which 
varies  depending  on  whether  there  is  an  even  or  an  odd  number  of  samples  (IV).  For 
N  odd, 

27 rif. 


(JV-l)/2 

Hdf{fs)odd  )  '  Oj  COS 

i=0 


fa 


For  N  even, 


H«U. 


N/2 

=  ^2  ai  cos 
i= 1 


(i  -  |)  fs 


fa 


(2.13a) 


(2.13b) 


For  these  MTFs,  a?-  is  the  digital  Liter  coefficient  and  fco  is  the  Liter  cut  oL  frequency 
in  cycles/mrad.  In  the  MatLab  FLIR92,  this  MTF  is  set  to  unity  for  all  frequencies 
in  both  the  horizontal  and  vertical  directions. 

k.  Display  MTF 

Although  many  recent  thermal  imagers  use  a  cathode  ray  tube  (CRT) 
display  to  communicate  collected  information  to  the  observer,  if  a  non-CRT  display  is 
being  used,  FLIR92  requires  a  user-specihed  MTF.  Otherwise,  the  model  calculates 
the  total  MTF  based  on  a  CRT  display.  In  FLIR92,  the  phosphor  spot  luminance 
intensity  is  assumed  to  have  a  Gaussian  distribution  with  a  o  value  either  specihed 
by  the  user  or  calculated  by  the  model  from  the  relationship 


a  = 


\ 


log(0.025) 
'  27 BpNr 


(2.14) 
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with  the  raster  frequency,  f^r  in  lines/mrad, 


(2.15) 


Here  the  model  uses  the  TV  limiting  resolution  criterion  that  the  MTF  equal  0.025  at 
the  raster  frequency.  Nr  is  the  number  of  lines  on  the  CRT  that  are  active.  (3  is  the 
vertical  FOV  in  mrad.  Using  these  quantities,  the  display  MTF  may  be  expressed  as 


(2.16) 


Since  this  result  is  only  valid  for  the  specified  resolution  criteria,  modifications  for 
other  resolution  criteria  may  be  made  by  multiplying  a  calculated  a  by  an  appropriate 
conversion  factor  and  then  entering  the  result  directly  into  the  model.  The  IREO 
Handbook  Vol.5  (1993)  again  points  out  that  a  follows  Equation  2.2b.  In  the  FLIR92 
MatLab  version  used  in  this  report,  a  was  hard  coded  to  0.01  mrad,  with  the  resultant 
MTF  shown  in  Figure  7. 
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Figure  7.  Display  MTF. 

The  standard  deviation  of  the  Gaussian  distribution  is  set  to  0.01  mrad. 
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1.  CCD  Charge  Transfer  Efficiency  MTF 

For  a  CCD  system,  Sequin  and  Thompsett  (1975)  explained  the  charge 
transfer  efficiency  MTF  as  due  to  the  non-unity  charge  transfer  efficiency  of  the 
system.  The  CCD  charge  transfer  efficiency  MTF  depends  on  the  number  of  gates 
in  the  transfer  from  the  detector  to  the  output  amplifier  (. N )  and  the  charge  transfer 
efficiency  at  each  gate  (e)  along  with  the  Nyquist  frequency,  or  the  sampling  frequency 
of  the  structure.  As  indicated  in  the  IREO  Handbook  Vol.4,  it  may  be  expressed  as 

HccdUs)  =  exp^-N(l  -  e)  1  -  cos  (^r^J  j  (2-17) 


where  rjv  is  again  the  Nyquist  frequency  in  cycles/mrad  (see  Section  II.A.l.d).  N 
may  be  calculated  in  the  model  by  assuming  that  the  system  has  a  interline  transfer 
scheme,  but  often  is  specified  by  the  user.  In  the  MatLab  program  being  used,  the 
CCD  Charge  Transfer  Efficiency  MTF  was  set  to  unity  for  all  frequencies,  in  both 
the  vertical  and  horizontal  directions. 

m.  Display  Sample  and  Hold  MTF 

Given  that  5S  is  the  sampling  aperture  in  mrad,  the  display  sample  and 


hold  MTF  has  the  form 


Hdsh{fs)  = 


sin  (7 rSsfs 


(2.18) 


Figure  8  shows  the  horizontal  display  sample  and  hold  MTF,  but  the  vertical  is  set 
in  the  MatLab  code  to  be  unity  for  all  frequencies. 

n.  Eye  MTF 

The  eye  MTF  encapsulates  the  influences  of  the  observer.  In  FLIR92, 
the  eye  MTF  is  based  on  work  by  Kornfeld  and  Lawson  (1971).  The  observer  may  be 
able  to  improve  and  optimize  the  system  display  by  adjusting  the  gain  and  level,  the 
viewing  distance,  and  the  like.  In  this  case,  the  eye  MTF  is  considered  non-limiting 
and  is  considered  a  constant  of  unity  since  no  degradation  of  the  spatial  frequency 
response  would  be  expected.  In  other  words, 


Heye{fs )  =  1-0. 


(2.19a) 
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On  the  other  hand,  the  observer  may  not  be  able  to  improve  or  adjust  the  display 
to  peak  performance.  If  this  is  the  case,  the  eye  limits  the  clarity  of  the  image,  so 
there  is  a  degradation  in  the  spatial  frequency  response.  The  MTF  to  describe  this 
degradation,  the  limiting  eye  MTF,  is  of  exponential  form, 

Heyeifs)  =  exp  '  (2.19b) 

Here,  M  is  the  system  magnification  and  V  is  the  spread  function  width  of  the  eye  in 
mrad.  F  depends  on  the  ambient  light  level  (Kornfeld  and  Lawson  1971).  Since  the 
observer  is  often  preoccupied  or  unable  to  adjust  the  display,  FLIR92  defaults  to  the 
limiting  eye  MTF.  In  the  MatLab  code,  V  is  hard  coded  to  2.196  mrad.  The  resultant 
MTF  is  shown  in  Figure  9. 

o.  System  Noise  Filter  MTF 

The  system  noise  filter  MTF  includes  only  the  detector  noise  and  is 
composed  of  all  the  noise  filter  MTFs.  These  MTFs  -  Hat,  Heh.p,  roll  up  of  the  temporal 
postfilter  MTF  (. HTPF ),  and  the  spatial  postfilter  roll  up  ( HSpf )  are  discussed  above. 


Sample  hold  MTF  H 


Figure  8.  Display  Sample  and  Hold  MTF  in  the  Horizontal  Direction. 
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The  system  noise  filter  MTF  is  not  included  in  the  total  system  MTF  ( Hsys )  since 
the  component  MTFs  are  already  included  in  the  total  MTF,  but  will  be  referenced 
in  future  sections.  It  may  be  expressed  as 

HNF(fs)  =  Hdi(fs)Hehp(fs)HTPF(fs)HSpF(fs )  (2.20) 

where  for  staring  systems,  HFpF(fs)  is  unity.  The  horizontal  system  noise  MTF  is 
shown  in  Figure  10. 

p.  Total  System  MTF 

The  total  system  MTF  is  the  product  of  all  the  component  MTFs  de¬ 
scribed  above.  For  the  generic  case  described  in  this  report,  the  total  horizontal  MTF 
is  shown  in  Figure  11.  The  downward  curve  of  the  MTF  illustrates  how  much  the 
image  is  degraded,  particularly  at  higher  spatial  frequencies,  corresponding  to  small 
targets. 


Display  brt  eye  response 


Figure  9.  Eye/Brain  Response  MTF. 
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2.  Sampling 

a.  Limits  to  Defined  MRTD 

In  FLIR92,  the  MRTD  is  only  defined  for  a  periodic  target  such  as  a 
4-bar  target.  The  target  must  have  a  7:1  aspect  ratio,  as  shown  in  Figure  1  and  the 
four  bars  should  be  fully  resolved  by  the  observer  in  order  to  specify  an  MRTD.  In 
thermal  imagers,  the  cutoff  frequency  for  the  observer  to  be  able  to  fully  resolve  the 
target  s  the  system’s  Nyquist  frequency.  By  definition,  then,  no  MRTD  can  be  given. 
This  limits  FLIR92  to  MRTD  prediction  below  the  Nyquist  limit. 

Since  many  targets  are  aperiodic,  this  artificial  method  of  limiting  the 
MRTD  prediction  may  be  pessimistic.  Information  relating  to  the  MRTD  may  be 
available  above  the  Nyquist  limit,  but  methods  for  obtaining  and  quantifying  this 
information  are  not  available  for  use  in  FLIR92. 


Horizontal  Noise  MTF 


Figure  10.  Horizontal  System  Noise  MTF. 
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b.  Sampling  Rate  and  Lattice  Spacing 

The  IREO  Handbook  Vol.4  (1993)  explains  that  imagers  sample  spa¬ 
tially  and  temporally.  Spatially,  the  thermal  imagers  sample  at  a  rate  of 

Rz  =  t-.  (2.21) 

This  is  half  the  Nyquist  frequency,  as  can  be  seen  by  comparing  Equations  2.5b  and 
2.21.  The  sample  lattice  spacing  is  then  the  inverse  of  Equation  2.21,  -2-.  FLIR92 
actually  calculates  the  effective  sampling  rate  and  effective  lattice  spacing  by  multi¬ 
plying  together  the  horizontal  and  vertical  components  (z  —  h  and  z  =  v)  and  then 
taking  the  square  root. 

c.  Optimizing  Scene  Phasing 

As  discussed  in  Section  II.A.l.d,  the  MRTD  measurement  may  be  de¬ 
graded  by  phasing  between  the  target  and  the  detector  grid.  There  is  an  MTF 
included  to  account  for  this  phase  difference,  particularly  when  dealing  with  staring 


Horizontal  System  MTF 


Figure  11.  Horizontal  System  MTF. 
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systems  where  the  MRTD  degradation  is  pronounced  when  there  is  a  misalignment 
between  respective  targets  and  detectors  under  consideration. 

3.  System  Noise 

FLIR92  characterizes  second-generation  thermal  imaging  systems  and  as  such, 
the  system  noise  in  MRTD  predictions  is  modeled  using  a  scaling  factor  that  multiplies 
the  random  spatio-temporal  noise  by  the  amount  of  excess  system  noise  (Kennedy 
1990).  To  simplify  the  model,  the  system  noise  is  reduced  to  components  that  add  in 
quadrature.  Appropriate  eye  spatial  and  temporal  integration  effects  are  considered 
for  each  noise  component  included. 

The  Night  Vision  and  Electronic  Sensors  Directorate  (NVESD)  of  the  US  Army 
in  1990  established  a  theoretical  framework  and  a  standard  laboratory  measurement 
procedure  to  better  characterize  the  noise  in  the  system  (Webb  et.al.  1991).  This 
method  of  noise  analysis  isolates  the  system  noise  into  eight  components.  These 
components  are  listed  in  Table  I  along  with  possible  sources  of  the  noise. 


Table  I.  3-D  Noise  Component  Descriptions  in  FLIR92 
Subscripts  of  noise  components  shown  in  table  indicate  dimension:  t  temporal,  v 
vertical  spatial,  h  horizontal.  Table  derived  from  Analyst’s  Reference  Guide  (1993). 


Noise 

Description 

Potential  Source 

(7tvh 

random  spatio-temporal  noise 

basic  detector  temporal  noise 

&tv 

temporal  row  noise 

line  processing,  j,  read-out 

(7th 

temporal  column  noise 

scan  effects 

(7vh 

random  spatial  noise 

pixel  processing,  detector-to-detector 
non-uniformity 

(Jy 

fixed  row  noise 

detector-to-detector  non-uniformity,  j 

<7h 

fixed  column  noise 

scan  effects,  detector-to-detector  non¬ 
uniformity 

(7t 

frame-to-frame  noise 

frame  processing 

s 

mean  of  all  components 

Here,  atvh  is  the  basic  detector  noise  which  is  often  characterized  by  the  NETD 
(see  Equation  2.22,  next  section).  The  NETD  is  a  sensitivity  parameter.  It  is  defined 
as  the  temperature  difference  between  the  target  and  the  background  that  will  produce 
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a  peak  signal  to  rms  noise  ratio  of  unity  at  the  output  of  a  reference  filter.  This  is  how 
the  NETD  gives  a  rough  estimate  of  how  noisy  a  sensor  is  compared  to  the  signals 
being  detected.  All  other  components  besides  atvh  must  be  added  in  quadrature 
in  order  to  completely  characterize  the  thermal  imaging  systems  noise  (Scott  et. 
al.  1993).  In  the  current  operational  FLIR92,  <jtvh  is  predicted  but  the  remaining 
components  are  measured  or  estimated.  FLIR92  includes  estimates  of  3-D  noise  for 
generic  scanning  and  staring  sensors  so  that  when  measurements  are  not  available, 
a  calculation  may  still  be  made.  More  detailed  descriptions  of  the  system  noise  are 
given  in  the  sections  following. 

a.  Defaults  for  3-D  Noise  Components 

In  a  given  system  where  3-D  noise  measurements  have  not  yet  been 
made  or  are  not  available,  FLIR92  provides  a  set  of  default  values.  These  default 
values  are  only  given  for  what  FLIR92  considers  significant  noise  components  and 
depend  on  the  predicted  random  spatial-temporal  noise,  atvh-  As  discussed  in  Section 
II. A. 3,  the  system  noise  sources  add  in  quadrature,  so  only  the  most  significant  noise 
components  are  defaulted  to  non-zero  components  since  the  others  would  scarcely 
affect  the  outcome  of  the  noise  calculation.  These  default  values  are  from  a  database 
of  system  noise  measurements  that  were  carried  out  by  NVESD  starting  in  April  1990. 
The  values  for  each  system  in  the  database  are  normalized  to  <jtvh  and  averaged  with 
other  systems  in  the  same  class  in  order  to  determine  which  are  the  dominant  noise 
sources  in  that  class  and  which  are  so  small  they  may  be  neglected. 

In  scanning  systems,  the  significant  noise  components  are  the  temporal 
row  and  the  fixed  row  noises,  atv  and  av.  Scanning  systems  typically  show  a  wide 
range  of  variation  in  noise  levels,  so  three  default  values  are  provided  for  both  atv  and 
<jv.  Table  II  shows  the  model’s  default  values  for  scanning  systems. 

In  staring  systems,  the  significant  noise  component  is  the  random  spa¬ 
tial  noise,  avh.  NVESD  found  a  single  default  value  for  avh  to  be  sufficiently  repre¬ 
sentative.  This  default  value  is  shown  in  Table  III. 
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Table  II.  3-D  Noise  Components  For  Scanning  Systems 
Subscripts  of  noise  components  shown  in  table  indicate  dimension:  t  temporal,  v 
vertical  spatial,  h  horizontal.  Table  derived  from  the  Analyst’s  Reference  Guide 
_ (1993). _ 


Noise  Term 

Low  Noise  Default 

Moderate  Noise  Default 

High  Noise  Default 

&tv 

0.25  atvh 

0.75  atvh 

1-0(7  tvh 

(Jv 

0.25atvh 

0.75  crtvh 

1-0(7  tvh 

Table  III.  3-D  Noise  Component  For  Staring  Systems 
Subscripts  of  noise  components  shown  in  table  indicate  dimension:  t  temporal,  v  ver¬ 
tical  spatial,  h  horizontal.  Table  derived  from  the  Analyst’s  Reference  Guide  (1993). 


Noise  Term 

Noise  Default 

&vh 

0A0atvh 

b.  Random  Spatio-temporal  Noise,  atvh 
Since  atvh  is  related  to  the  actual  bandwidth  of  the  system,  it  may  be 
measured  directly  at  the  output  port  prior  to  display.  However,  as  outlined  in  the 
IREO  Handbook  Vol.4,  it  may  otherwise  be  determined  through  the  relationship  to 


the  NETD, 


<?tvh  =  NET D  x 


A  fp 


A/. 


(2.22) 


N 


where  A is  the  equivalent  noise  bandwidth  for  the  NETD.  In  this  case,  A fp  is 
the  actual  noise  bandwidth  which  is  associated  with  the  system  electronics  before 
display.  In  order  to  determine  the  equivalent  noise  bandwidth,  given  that  S(v)  is 
the  normalized  detector  noise  power  spectrum  and  Href(is )  is  the  standard  NETD 
reference  filter,  the  relationship  used  is 


A/at  =  J  S{v)[Href(v)]2  dv.  (2.23) 

o 

The  actual  noise  bandwidth  differs  depending  on  whether  a  scanning  or  staring  system 
is  under  consideration.  When  considering  scanning  systems,  the  noise  bandwidth  at 
the  system  output  port  is  of  similar  form  to  the  equivalent  noise  bandwidth, 


A/p  =  J  S(n)[HTPF(o)}2dn.  (2.24) 

o 
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For  staring  systems,  the  actual  noise  bandwidth  is  given  by 


(2.25) 


and  t,  is  the  FPA  integration  time. 

For  all  classes  of  systems,  staring  or  scanning,  the  general  form  of  the 
random  spatio-temporal  noise  is  given  in  Equation  2.26, 


®tvh 


.  ^2 


(2.26) 


7rW^/£>*(A,300)^(A)dA 

Ai 


Here,  fno  is  the  optical  f-number,  A fp  is  the  system  noise  bandwidth  as  defined 
in  Equation  2.24  or  2.25,  To  is  the  optical  transmittance,  and  A ^  is  the  detector 
area  in  cm2.  The  partial  derivative  is  the  thermal  derivative  of  Planck’s  Law  in 
W/crir/sr/pm.  The  detector  noise- limited  spectral  detectivity  is  D*(A,300)  in  cm- 
Hz^ /W  and  includes  only  noise  components  from  the  temporal  noise  sources  since 
the  spatial  noise  source  contributions  are  included  in  the  system  noise  correction 
functions  described  in  the  following  section,  I1.A.3.C.  Also,  FL1R92  does  not  make 
any  adjustments  for  changes  to  detector  responsivity  due  to  cold  shielding,  so  any 
tweaking  there  must  be  off-line. 

c.  Noise  Correction  Functions 

Since  each  3-D  noise  component  listed  in  Table  1  is  assumed  to  be 
statistically  independent,  the  total  system  noise,  filtered,  may  be  written  as  the  root 
sum  square  of  the  noise  components: 


«(/.)  =  °lhE,E„SmKU,)  +  °lkEv,U,)EhSf,) 

+  c2thEtEh,(f,)  +  aiE,Ev,if.)  +  crlE„,(f.) 

+  °2hEh,(f,).  (2.27) 


Here,  Et,  Eflz(fs)  and  EVz(fs)  represent  the  eye  or  brain  temporal  and  spatial  integra¬ 
tion  associated  with  the  noise  components.  The  frame-to-frame  term  of  FA  has  been 


dropped  because  the  frame- to- frame  noise  at  is  almost  always  so  small  comparatively 
that  it  is  negligible.  The  orientation  of  the  MRTD  target  under  consideration  is  in¬ 
dicated  by  the  subscript  z.  In  this  case,  the  temporal  integration  may  be  expressed 
in  terms  of  Fr,  the  system  frame  rate  in  Hz,  te,  the  eye  integration  time  in  sec,  and 
at,  the  temporal  sample  correlation  length. 

Ft  =  "yr~  •  (2.28) 

r  rTe 

In  FLIR92,  at  is  assumed  to  be  unity  in  MRTD  calculations.  The  spatial  integra¬ 
tions  are  simplified  in  the  FLIR92  model  and  may  be  expressed  in  terms  of  sample 
correlation  factors  ah  and  av,  horizontal  and  vertical  sampling  rates  Rh  and  Rv  in 
samples/mrad,  and  spatial  integration  limits  Lhz(fs)  and  LVz(fs)  in  mrad-1.  These 
spatial  integration  limits  are  approximately  the  width  and  height  of  the  MRTD  bar 
target.  For  staring  systems,  the  samples  in  each  direction  are  assumed  independent 
so  that  ah  and  av  are  unity.  For  scanning  systems,  samples  in  the  scan  direction 
may  not  be  assumed  independent  due  to  the  motion  of  the  scanner,  so  ah  may  be 
greater  than  unity.  Similarly,  samples  taken  perpendicular  to  the  scan  direction  may 
be  assumed  independent  since  the  motion  of  the  scanner  is  cross-directional,  so  av 
is  unity.  Although  FLIR92  uses  this  simplified  form,  the  exact  expressions  for  the 
horizontal  and  vertical  eye/brain  spatial  integration  are  given  in  Appendix  B. 

In  order  to  determine  the  noise  terms  for  the  horizontal  and  vertical 


MRTD  calculations,  it  is  assumed  that  only  noise  components  in  the  direction  of  the 
MRTD  target  degrade  the  MRTD.  For  the  MDTD,  since  target  orientation  does  not 
affect  the  calculation,  the  noise  correction  function  is  independent  of  direction.  Given 


the  noise  correction  functions  in  the  horizontal, 


and  in  the  vertical, 


(2.29) 


(2.30) 
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the  horizontal  and  vertical  MRTD  noise  terms  may  be  expressed.  In  the  horizontal 
direction,  the  noise  term  is  given  as 


Vh(fs)  =  Vtvhkh{fs)jEtEVh(fs)Ehh(fs 


(2.31) 


and  in  the  vertical, 


wv{fs)  =  Vtvhkv(fs)jEtEVv(fs)Ehv(fs). 


(2.32) 


For  the  MDTD,  the  noise  correction  function  has  the  form 


n(/s)  —  1  + 


2  2 
® Vh  TTI—  1  i  ^tl 

~2 —Et  +  -f 


f-  [^(/.)]_1  +  [EtEv(fs)Y 


[Eh(f.)]-1  +  -%-[EtEh(f.)]- 


4.  Calculating  the  Predicted  MRTD 

The  basic  form  for  the  MRTD  calculated  in  FLIR92  is  of  the  form 


(2.33) 


MRTDz(fs)  = 


\SNRTHatvhkz(fs 

H  SYS,  if  s') 


/EtEhz(fs)EVz(fs 


(2.34) 


where  z  again  indicates  the  orientation  of  the  MRTD  target  in  either  the  horizontal 
or  the  vertical  direction.  The  denominator  HsYSz(fs)  is  the  total  system  MTF,  or 
the  product  of  all  component  MTFs.  kz(fs )  is  the  3-D  noise  correction  function,  as 
given  in  Equations  2.29  and  2.30.  The  threshold  SNR  required  to  resolve  the  MRTD 
target  is  SNRTH  and  the  eye/brain  integrals  are  as  described  in  Section  II. A. c  and 
Appendix  B.  The  IREO  Handbook  Vol.5  further  outlines  this  method  of  determining 
the  predicted  MRTD.  More  details  about  the  included  variables  in  this  general  form 
of  the  MRTD  calculation  are  given  below.  The  system  MRTD  for  the  given  input 
parameters  is  shown  in  Figure  12. 

a.  Psychophysical  Constants  SNRth  and  te 
In  FLIR92,  the  threshold  SNR  to  resolve  the  MRTD  target  and  the  eye 
integration  time,  SNRth  and  te,  may  be  adjusted  to  tweak  an  MRTD  prediction  to 
better  match  a  set  of  measurements.  NVESD  does  have  recommended  settings,  as 
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given  in  Table  IV.  Biberman  (1973)  showed  using  psychophysical  data  that  SNRTH 
is  a  function  of  the  target  spatial  frequency,  but  the  NVESD  recommended  value  of 
2.5  is  a  representative  average  for  optimal  observing  conditions. 


Table  IV.  NVESD  Recommended  Settings  for  Psychophysical  Constants 
Settings  in  FLIR92  may  be  adjusted  from  NVESD  recommendations. 


Psychophysical  Constant 

NVESD  Recommended  Setting 

Units 

SNRth 

2.5 

- 

te 

0.1 

s 

Te  depends  on  the  background  luminance.  Luminances  corresponding 
to  a  0.1  s  eye  integration  time  show  agreement  with  display  luminances  that  NVESD 
measured  in  perception  experiments  in  1988.  These  experiments  were  conducted 
under  conditions  similar  to  those  used  in  MRTD  measurements  such  as  a  darkened 
room  and  optimal  viewing.  Observers  set  the  display  luminance  to  an  average  of 
0.15  milli-Lamberts.  Higher  ambient  light  conditions  typically  correspond  to  greater 


Test  horizontal  and  vertical  MRT 


Figure  12.  System  MRTD  showing  the  Horizontal  and  Vertical  Results. 
The  2-D  system  MRTD  is  calculated  by  taking  the  geometric  average  of  the 
horizontal  and  vertical  results,  so  it  would  lie  between  the  two. 
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display  luminances  and  a  faster  integration  time  may  be  appropriate  in  these  cases, 
such  as  in  the  systems  used  in  the  held. 

b.  2-D  MRTD  Calculation 

FLIR92  can  calculate  a  2-D  MRTD  by  taking  the  geometric  mean  of 
the  horizontal  and  vertical  MRTDs.  In  taking  the  geometric  mean,  each  component 
is  weighted  equally  and  the  mean  is  with  respect  to  the  spatial  frequency  axis.  This 
causes  the  2-D  MRTD  to  asymptote  at  the  mean  value  of  the  vertical  and  horizontal 
cut-off  spatial  frequencies.  These  cut-off  spatial  frequencies  are  determined  by  either 
the  Nyquist  limits  or  MTF  roll  off. 

5.  Calculating  the  Predicted  MDTD 

The  basic  form  for  the  MDTD  calculated  in  FLIR92  is 


MDTD(f3)  =  jEtEkU,)EM.)  (2.35) 

where  the  MDTD  is  independent  of  target  orientation,  so  /cmdt(/.s)  is  independent 
of  orientation  and  is  given  by  Equation  2.33.  Qh(fs),  Qv(fs),  and  Ev(fs )  are  defined 
below  and  are  equivalent  for  both  scanning  and  staring  systems.  Along  the  horizontal 
direction,  the  eye/brain  spatial  integration  differs  for  staring  and  scanning  systems, 
so  each  is  stated  below.  At  is  the  target  area  in  mrad2  and  is  related  to  the  spatial 


frequency,  fs  by 


At  =  — 


(2.36) 


Equation  2.36  assumes  an  isometric  target,  but  is  a  poor  approximation  for  many 
operation  targets  that  are  not  so  uniform.  All  other  variables  are  defined  as  for  the 
MRTD  calculation. 

The  Qh(fs)  and  Qv(fs )  integrals  both  include  the  total  system  MTF,  Hsys- 
If  z  represents  the  orientation,  either  horizontal  or  vertical,  then  the  Qz{fs)  integral 


Qz(fs)  —  /  [ hSysz(v)Y 


(2.37) 
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Ev(fs) ,  the  vertical  eye/brain  spatial  integration,  is  defined  the  same  way  in  both  the 
staring  and  the  scanning  systems  and  is  of  the  form 


Ev{fs)  —  —  l  [Hsysv (y)}2 [HNfv (v)Y 


sm 


(f) 


7 TV 
fs 


dv. 


(2.38) 


Eh(fs)  is  defined  differently  for  staring  and  scanning  systems.  The  coefficient  in  front 
of  the  integral  changes.  For  scanning  systems  it  is  of  the  form, 

(tT2 


Eh{fs)  =  ^r  J  S{v)[HSYSh(v)Y\HNFh(v)Y 


sm 


■KV 

fs 


dv. 


(2.39) 


For  staring  systems,  it  is  similar  to  the  vertical  direction  given  in  Equation  2.38, 


r  O' 

Eh(fs)  =  ~  J  {HSYSh(v)]2[HNFh(v)Y 


sm 


(f) 


7 TV 
fs 


dv. 


(2.40) 


6.  MRTD  and  MDTD  Temperature  Dependence 

Since  background  temperature  affects  the  random  spatial-temporal  noise  fac¬ 
tor,  atvh ,  both  the  MRTD  and  MDTD  calculations  are  influenced  by  the  background 
temperature.  In  FLIR92,  the  default  background  temperature  is  300  K,  but  the  model 
can  scale  atvh  to  any  background  temperature  between  240  K  and  330  K.  Generally, 
though,  predictions  made  with  a  background  temperature  of  300  K  are  acceptable  for 
most  FLIR92  applications.  In  fact,  atvh  is  the  only  3-D  noise  component  that  FLIR92 
allows  to  be  scaled  to  variable  background  temperatures.  Temperature  scaling  only 
takes  into  account  the  background  temperature;  there  is  no  scaling  to  account  for  the 
temperature  dependence  of  transmittance  in  optical  materials  or  noise  from  photon 
emissions  by  variably  temperatured  sensor  components.  For  details  of  the  tempera¬ 
ture  scaling  of  atvh,  see  Appendix  C. 


7.  Johnson  Criteria  and  FLIR92 

The  Johnson  criteria  uses  the  resolving  power  of  an  imager  under  consideration 
to  provide  an  estimate  for  how  good  the  imager  will  be  at  target  acquisition.  The 
resolving  power  of  the  imager  in  cycles/mrad  is  multiplied  by  the  target  size  in  mrad  to 
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get  the  number  of  cycles  on  target.  This  number  of  cycles  on  target  is  then  compared 
to  a  table  of  the  number  of  cycles  on  target  needed  to  detect,  recognize,  identify,  or 
otherwise  perform  target  acquisition  tasks. 

B.  NVTHERMIP 

Prior  to  the  introduction  of  NVThermlP,  NVTherm  was  developed  to  rectify 
known  errors  in  FL1R92,  particularly  as  regards  staring  arrays.  The  FL1R92  model 
did  not  take  into  account  the  contrast  limitations  of  the  observer’s  eye.  Also,  FLIR92 
limitations  on  detector  size,  spacing,  and  fill  factor  caused  under-sampled  imagery. 
The  new  model  attempted  to  compensate  for  the  sampling  artifacts  by  imposing  a 
cutoff  at  half  the  sample  rate  of  the  imager,  but  when  this  result  was  used  with  the 
Johnson  metric,  resultant  predictions  were  pessimistic.  Although  the  basic  MRTD 
theory  used  in  FLIR92  went  into  NVTherm  and  NVThermlP  largely  unchanged,  NV- 
ThermlP  reworked  the  eye  MTF  contribution  to  improve  the  pessimistic  results  from 
FL1R92.  NVThermlP  predicts  the  CTF  of  thermal  imaging  systems,  but  can  also 
predict  a  laboratory  measurement  of  the  the  MRTD.  NVThermlP  also  has  completely 
redone  the  eye  MTF.  In  this  report,  the  CTF  is  used  for  comparison  to  the  MRTD 
so  the  new  eye  MTF  theory  will  be  included  in  the  comparison. 

1.  MTFs  Included  in  NVThermlP 

The  MTFs  included  in  NVThermlP  are  predominantly  the  same  as  those  in¬ 
cluded  in  FLIR92.  As  in  FLIR92,  in  most  cases  a  user-defined  MTF  may  be  entered 
and  is  encouraged  whenever  possible,  for  more  realistic  predictions.  All  the  follow¬ 
ing  MTFs  are  derived  from  the  NVThermlP  User  Manual  (NVESD  2005)  unless 
otherwise  stated. 

a.  Optical  MTF 

The  optical  MTF  still  includes  the  diffraction- limited  MTF  and  geo¬ 
metric  blur  MTF,  both  of  the  same  forms  as  Equations  2.1  and  2.2a.  NVThermlP 
assumes  the  optical  MTF  is  the  same  in  the  horizontal  and  vertical  directions. 
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b.  Detector  Spatial  MTF 

NVThermlP  calculates  the  detector  spatial  MTF  using  Equation  2.3b 
and  does  not  distinguish  between  SPRITE  and  non-SPRITE  systems.  A  user-specified 
MTF  would  be  required  to  clarify.  Hds{fs)  is  not  assumed  to  be  symmetric  in  the 
horizontal  and  vertical  directions. 

c.  Focal  Plane  Array  Integration  MTF 

The  focal  plane  array  integration  MTF  is  not  included  in  NVThermlP, 
but  could  be  entered  by  the  user  as  one  of  the  custom  MTFs  if  so  desired. 

d.  Sample-scene  Phase  MTF 

The  sample-scene  phase  MTF  is  not  included  in  NVThermlP,  but  could 
be  entered  by  the  user  as  one  of  the  custom  MTFs  if  so  desired. 

e.  Image  Motion  MTF 

In  NVThermlP,  the  only  image  motion  MTF  included  is  the  random 
motion  MTF  that  describes  the  vibration  of  the  sensor  platform.  As  in  FLIR92 
Equation  2.7,  Hmr(fs)  is  assumed  to  be  Gaussian.  It  is  not  assumed  symmetrical  in 
both  directions.  User-defined  MTFs  for  linear  and  sinusoidal  motion  could  be  entered, 
or  a  total  motion  MTF  defined. 

/.  Detector  Temporal  MTF 

The  detector  temporal  MTF  is  not  included  in  NVThermlP,  but  could 
be  entered  by  the  user  as  one  of  the  custom  MTFs  if  so  desired. 

g.  Electronics  Frequency  Responses 

NVThermlP  calculates  a  low  pass  filter  MTF  for  scanning  sensors,  only 
in  the  horizontal  direction.  The  MTF  is  the  same  as  Equation  2.10b.  No  high  pass 
filter  MTF  is  outputted. 

h.  Boosting  MTF 

The  boosting  MTF  is  not  included  in  NVThermlP,  but  could  be  entered 
by  the  user  as  one  of  the  custom  MTFs  if  so  desired. 
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i.  Electro- Optical  Multiplexer  MTF 

In  NVThermlP,  the  electro-optical  multiplexer  MTF  is  calculated  by 
Equation  2.12  unless  the  user  specifies  an  MTF.  It  varies  in  the  horizontal  and  vertical 
directions  due  to  the  possible  differences  in  the  LED  size. 

j.  Digital  Filter  MTF 

The  digital  filter  MTFs  in  NVThermlP  are  similar  to  those  in  FLIR92, 
Equations  2.13a  and  2.13b,  but  have  a  factor  to  allow  for  dither.  The  dither  factor, 
D,  is  set  to  1  when  the  user  chooses  not  to  include  dither.  The  modified  equations 
are  given  below.  For  N  odd  the  MTF  is 

(AT-l)/2 


Hdf(fs)odd  ^  ) 


ai  cos 


i= 0 


and  for  N  even, 


Had!, 


N/2 

=  ^2  ai  cos 

i= 1 


27 rifs 
Dfco 

^  0 -  0  f* 

Dfco 


(2.41a) 


(2.41b) 


k.  Display  MTF 

Though  the  display  MTF  in  NVThermlP  assumes  a  CRT  display  for 
presenting  information  to  the  observer,  it  includes  a  factor  Ezoom  that  is  not  in 
FLIR92,  Equation  2.16.  Ezoom,  Ezoom,  is  related  to  the  magnification  in  the  system 
parameters  and  this  magnification  is 

0  arctan 


M  = 


eye 


FOVq, 


=  2- 


\2DVd  J 


FOV„ 


(2.42) 


where  Dh  is  the  display  height  in  cm,  Dvd  is  the  display  viewing  distance  in  cm, 
and  FOVv  is  the  vertical  sensor  FOVvin  mrad.  When  Ezoom  is  used,  it  is  assumed 
that  only  part  of  the  full  FOV  will  be  seen  on  the  display  screen.  Ezoom  examples 
are  given  in  Table  V.  Values  for  M  may  vary  from  l /6th  to  200,  but  typically  are 
between  0.5  and  20. 

Equation  2.16  is  modified  to  include  the  Ezoom  factor  by  modifying  a: 


a  = 


SwFOVz 

F zoomD  H 


(2.43) 
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Table  V.  Ezoom  Examples 

Portions  of  total  FOV  area  shown  on  display  screen  when  different  Ezoom  values  are 
used.  Values  obtained  from  the  NVThermlP  User’s  Manual  (2005). 


Ezoom 

M  Factor 

Vert  FOV  Seen 

Horiz  FOV  Seen 

Tot  FOV  Area  Seen 

none 

1 

1 

1 

1 

single 

2 

1/2 

1/2 

1/4 

double 

4 

1/4 

1/4 

1/16 

where  Sw  is  the  physical  width  of  the  Gaussian  display  spot  in  cm  and  FOVz  is  the 
FOV  in  either  the  horizontal  (z  =  h )  or  the  vertical  (z  =  v)  direction. 

Also,  unlike  FL1R92,  NVThermlP  includes  MTFs  for  other  possible 
display  types  besides  CRT,  but  these  are  not  discussed  here  since  this  report  focuses 
on  comparing  the  two  models  using  the  same  inputs. 

l.  CCD  Charge  Transfer  Efficiency  MTF 

The  CCD  Charge  Transfer  Efficiency  MTF  is  not  included  in  NVThe¬ 
rmlP,  but  could  be  entered  by  the  user  as  one  of  the  custom  MTFs  if  so  desired. 

m.  Display  Sample  and  Hold  MTF 

The  display  sample  and  hold  MTF  of  Equation  2.18  is  again  the  same 
in  NVThermlP,  but  is  only  calculated  for  scanning  systems.  Also,  it  is  only  applied 
when  calculating  the  horizontal  MTF. 

n.  Eye  MTF 

NVThermlP  handles  the  eye  MTF  very  differently  from  FLIR92.  NV- 
ThcrmlP  considers  the  human  eye  point  spread  function  as  a  combination  of  three 
factors:  the  optics,  the  retina  and  the  tremor.  This  leads  to  an  MTF  that  may  be 
expressed  as  the  product  of  the  component  MTFs  from  the  eye  optics  He0)  the  retina 
Hret,  and  the  tremor  Htrem: 


Heye(fa )  =  Heo(fs)Hret(fs)Htrern(fs).  (2.44) 

Overington  (1976)  outlines  the  theory  of  the  human  eye  MTF  and  identifies  the  three 
component  MTFs  above.  Based  on  Overington’s  work,  the  general  forms  for  the  eye 
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MTF  are  found,  but  in  In  NVThermlP,  these  forms  for  the  eye  optical,  retinal,  and 
tremor  MTF  are  all  empirical  (NVESD  2001).  The  eye  optics  MTF  is 


=  exp 


/ 

Ja 

io\ 

_ 

43  6g  V  MW 

V 

fo 

J 

(2.45) 


where  Msys  is  the  imaging  system  magnification,  and  fQ  and  iQ  are  defined  below, 


f0  =  exp  (3.663  -  0.0216DJ  log10(£)p)) 


=  0.7155  + 


0.277 


(2.46a) 

(2.46b) 


The  variable  Dp  is  the  pupil  diameter  in  mm  and  is  valid  if  one  eye  is  used.  If  two 
eyes  are  used,  then  NVThermlP  reduces  the  pupil  diameter  by  0.5mm.  Dp  is  defined 
as 


Dp  =  —9.011  +  13.23  exp 


log  10  (A 


(2.46c) 


21.082 

\  / 

where  Jl  is  the  number  of  foot-Lamberts  at  the  eye  from  the  display  and  is  /l  =  ylfg- 
Ld  is  the  display  luminance  in  milli-Lamberts.  The  retina  MTF  is  defined  as 

l  1.210 


Hret{fs)  =  exp  -0.375 


fs 


M, 


sys 


(2.47) 


The  tremor  MTF,  or  the  MTF  of  the  eye  due  to  tremor,  is 


HtremUs )  =  exp  -0.4441 


fs 


M, 


sys 


(2.48) 


Clearly,  then,  the  MTFs  depend  on  the  pupil  diameter,  which  in  turn  depends  on  the 
light  level. 

Figure  13  compares  the  eye  MTFs  in  FLIR92  and  NVThermlP.  Clearly, 
NVThermlP  has  a  more  optimistic  estimate  of  the  image  degradation  due  to  the 
observer’s  eye,  which  reduces  the  known  error  in  the  FLIR.92  modeling  of  the  human 
eye. 
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o.  System  Noise  Filter  MTF 

The  system  noise  filter  MTF  is  merely  the  roll-up  of  all  the  system 
noise.  In  NVThermlP,  there  is  a  further  factor  for  system  noise  that  is  requested  as  a 
user-input.  NVThermlP  considers  that  each  system  has  a  different  fixed  pattern  noise 
associated  with  variations  in  detector  gain  and  level  offset.  There  are  three  options 
in  the  NVThermlP  input:  None,  Noise  Factor,  and  3-D  Noise.  If  None  is  selected, 
then  the  system  is  considered  ideal  with  no  variations  in  gain  and  level  offset  among 
the  detectors.  For  Noise  Factor ,  experimentally-determined  independent  horizontal 
and  vertical  factors  are  multiplied  by  the  horizontal  and  vertical  noise  CTFs.  The 
3-D  Noise  is  discussed  in  a  later  section  about  system  noise. 

p.  Interpolation  MTF 

NVThermlP  includes  an  MTF  for  interpolation,  or  the  process  of  in¬ 
creasing  the  image  size  by  inserting  ‘filler’  pixels  between  the  original  pixels.  Interpo¬ 
lation  may  be  done  either  vertically  or  horizontally,  and  if  interpolation  is  selected  for 


Compare  Eye  MTF  from  FLIR92,  NVThermlP 


Figure  13.  Comparing  the  Eye  MTF  in  FLIR92  and  NVThermlP  in  the  Horizontal 
Direction. 
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both  directions,  only  one  is  applied  at  a  time.  The  results  are  then  combined.  NV- 
ThcrmlP  offers  three  different  methods  for  interpolation,  Pixel  Replication ,  Bilinear 
Interpolation,  and  Bicubic  Interpolation.  In  this  report,  since  FLIR92  does  not  offer 
interpolation,  no  interpolation  is  selected. 

q.  Optical  Turbulence  MTF 

Unlike  FLIR92,  operational  NVThermlP  already  includes  an  MTF  for 
optical  turbulence.  This  turbulence  MTF  assumes  that  C‘ZI  is  constant  along  the 
pathlength  and  so  is  an  average  MTF.  is  a  parameter  that  describes  the  optical 
turbulence  in  the  atmosphere  and  more  thoroughly  defined  in  Chapter  III.  The 
formulation  is  very  similar  to  the  Goodman  results  for  constant  given  in  Equation 
A. 98  in  Appendix  A.  The  equation  is  given  below, 

Hat(fs)  =  exp  ^—57 Aafi C^\~^ z  ^1  -  \  j  j  (2.49) 

where  a  is  a  constant  defined  as  |,  fs  is  here  defined  in  cycles/rad,  C\  is  still  in  m~i, 
and  A,  z,  and  D  are  all  in  m. 


2.  Sampling 

a.  Sample  Spacing 

For  an  imaging  system,  the  sample  spacing  quantifies  the  limits  due  to 
sampling.  For  staring  systems,  the  sample  spacing  can  be  calculated  in  the  horizontal 
and  vertical  directions  by 

&ssz  =  (2.50a) 

where  z  indicates  the  target  orientation,  either  horizontal  or  vertical.  Nz  is  the  number 
of  detectors  in  the  z  direction.  assz  is  in  mrad.  For  a  scanning  system,  the  vertical 
sample  spacing  is  calculated  as  per  Equation  2.50a,  but  in  the  horizontal  direction 
the  sample  spacing  is  different.  There  is  no  sample  spacing  in  a  continuously  scanning 
system,  so  assh  =  0.  A  scanning  system  that  samples  has  a  sample  spacing  requires 


a  user-input  of  NHIFOv,  or  samples  per  HIFOV.  Then  the  sample  spacing  is 


assH 


DASh 

Nhifov 


(2.50b) 
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DASh  is  the  detector  angular  subtense  (DAS)  in  the  horizontal  direction,  in  other 
words,  the  detector  width  divided  by  the  focal  length  of  the  collecting  optics. 

b.  Sample  Frequency 

As  explained  in  Section  lI.A.2.b,  the  sampling  frequency  is  just  the 
inverse  of  the  sample  spacing,  Equations  2.50a  or  2.50b.  Thus,  the  sampling  frequency 

is 


f samp 


aSSz 


(2.51) 


Since  the  horizontal  and  vertical  sample  spacing  may  differ,  the  horizontal  and  vertical 
fsamp  may  also  differ.  Again  the  half  sample  frequency  is  the  Nyquist  frequency. 


3.  System  Noise 

As  mentioned  in  Section  II.B.l.o,  the  3-D  system  noise  is  handled  in  a  similar 
manner  to  FL1R92.  The  3-D  system  noise  components  are  defined  in  Table  VI, 
essentially  the  same  as  Table  I. 


Table  VI.  3-D  Noise  Component  Descriptions  in  NVThermlP 
Subscripts  of  noise  components  shown  in  table  indicate  dimension:  t  temporal,  v  ver¬ 
tical  spatial,  h  horizontal.  Table  derived  from  NVThermlP  User’s  Manual  (NVESD 
2005). 


Noise 

Description 

Potential  Source 

&tvh 

random  spatio-temporal  noise 

Basic  Detector  Temporal  Noise 

&tv 

temporal  row  noise,  line  bounce 

Line  Processing,  j,  read-out 

(7  th 

temporal  column  noise,  column 
bounce 

Scan  Effects 

(7yh 

random  spatial  noise,  bi¬ 
directional  fixed  pattern  noise 

Pixel  Processing,  Detector-to-Detector 
Non-Uniformity  j 

(Jy 

fixed  row  noise,  line-to-line  non¬ 
uniformity 

Detector-to-Detector  Non-Uniformity 

<7h 

fixed  column  noise,  colunm-to- 
column  non-  uniformity 

Scan  Effects,  Detector-to-Detector 
Non- LI  niformity 

<7t 

frame-to-frame  noise,  frame 
bounce 

frame  processing 

s 

mean  of  all  noise  components 
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As  in  FLIR92,  these  eight  noise  parameters  are  derived  using  directional  av¬ 
erages.  atvh  is  predicted  and  all  other  components  are  estimated  based  on  historical 
databases  of  measurements.  In  Table  VI,  the  subscript  that  is  missing  indicates  which 
directions  were  averaged  so,  for  example,  crtv  was  averaged  in  the  horizontal  and  av 
in  both  the  temporal  and  horizontal.  Again  similar  to  FLIR92,  only  certain  noise  pa¬ 
rameters  are  considered  important  in  scanning  and  staring  systems.  Table  VII  shows 
the  key  parameters  for  scanning  systems,  where  the  noise  term  is  normalized  using 
the  random  spatio-temporal  noise  crtvh. 

Table  VII.  Noise  Values  for  Scanning  Systems  in  NVThermlP 
Subscripts  of  noise  components  shown  in  table  indicate  dimension:  t  temporal,  v  ver¬ 
tical  spatial,  h  horizontal.  Table  derived  from  NVThermlP  User’s  Manual  (NVESD 
2005). 


Noise  Term 

Low 

Moderate 

High 

@vh/  &tvh 

0 

0 

0 

/ &tvh 

0.25 

0.75 

1 

@tvh 

0 

0 

0 

Table  VIII  shows  the  key  parameters  for  staring  systems,  where  the  noise  term 
is  again  normalized  with  crtvh-  Table  VIII  compares  to  Table  III. 


Table  VIII.  Noise  Values  for  Staring  Systems  in  NVThermlP 
Subscripts  of  noise  components  shown  in  table  indicate  dimension:  t  temporal,  v  ver¬ 
tical  spatial,  h  horizontal.  Table  derived  from  NVThermlP  User’s  Manual  (NVESD 
2005). 


Noise  Term 

Low 

Moderate 

High 

@vh/ &tvh 

0.2 

0.5 

1  -  2 

/  ®tvh 

0.2 

0.5 

1  -  2 

@h/ &tvh 

0.2 

0.5 

1  -  2 
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a.  Random  Spatio-temporal  Noise,  atvh 

The  random  spatio-temporal  noise  is  calculated  similarly  in  NVThe- 
rmlP  as  in  FLIR92,  but  includes  a  factor  to  account  for  the  number  of  detectors. 
Again,  atvh  is  calculated  assuming  a  background  temperature  of  300K  and  the  rela¬ 
tionship, 


&tvh  — 

Ai 


(2.52) 


where  N  is  the  number  of  detectors  and  all  other  variables  are  defined  as  for  Equation 
2.26.  Note  that  Equation  2.52  differs  from  Equation  2.26  by  a  factor  of  n  in  the 
denominator. 


4.  Calculating  the  Predicted  MRTD 

NVThermlP  outputs  three  laboratory  measurements  of  the  MRTD:  high  gain, 
low  gain,  and  user  input.  The  basic  equation  used  in  all  cases  is 


MRTDz(fs) 


2  SCNTMPCTFeye(fs)  a2Y2detQHz(fs)QVz(fs ) 

(Abarifs)  -  Aspc(fs))SL  \  SCNLp 


(2.53) 


Here,  SCNTMP  is  the  scene  contrast  temperature  that  generates  the  average  display 
luminance  in  K,  CT Feye(fs)  is  the  naked  eye  CTF,  Abar  is  the  area  of  the  bar  in  the 
target  in  cm2,  Aspc  is  the  area  of  the  space  between  bars  in  the  target  in  cm2,  Sp  is 
the  unitless  normalized  laboratory  detector  responsivity,  Y^et  is  the  noise  variance  in 
K-mrad-s^ ,  and  QHZ  and  QVZ  are  the  horizontal  and  vertical  noise  bandwidth  for  the 
horizontal  system  CTF  and  the  vertical  system  CTF.  The  laboratory  conditions  can 
be  strictly  controlled,  which  is  important  because  in  this  calculation  when  Abar(fs )  = 
Aspc{fs),  the  MRTD  will  be  undefined. 

The  laboratory  MRTD  in  NVThermlP  is  difficult  to  compare  directly  to  the 
MRTD  in  FLIR92  for  many  reasons.  Two  of  particular  concern  in  this  study  were 
the  high,  low,  and  user  gain  of  NVThermlP  -  there  is  no  equivalent  in  FL1R92  -  and 
the  fact  that  NVThermlP  calculates  the  MRTD  for  laboratory  conditions  which  are 
highly  controlled  and  not  equivalent  to  the  conditions  assumed  in  FL1R92.  Since 
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the  laboratory  MRTD  outputted  by  NVThermlP  is  not  directly  comparable  to  the 
MRTD  from  FLIR92,  a  work-around  was  proposed  in  conversations  with  Jon  Hixson 
(Army  Research  Labs,  Fort  Belvoir,  VA).  The  CTF  outputted  from  NVThermlP  can 
be  manipulated  to  obtain  an  MRTD  comparable  to  FLIR92  by  multiplying  2  x  the 
SCT.  NVThermlP  produces  a  system  CTF  in  the  horizontal  and  vertical  directions 
following  Equations  2.54a  and  2.54b: 


CTF„,Jf,h) 

CTFVn,(f„) 


\ 


{CTFH,Jf,„)P  + 


CTF„  (/. 


Shi 


SCNi 


tmp 


\ 


{CTFv(f„)]i  + 


CTFV  .  (/, 


SCNi 


tmp 


(2.54a) 


(2.54b) 


where  the  eye  CTFs  and  the  noise  CTFs  in  the  horizontal  and  vertical  directions  are 
defined  as 


erFi/noisc(/sj 

CTFVnoi,'(fsv) 

CTFH,Jf„) 

CTFV,JU 


(  CTFeye 

(  fsh  \  N 

V Smag ) 

(2.56a) 

\MdlspMTFHsys(fSh) f 

(  CTFeye  ( 

'  fsv  )  \ 

\] a2^detQHverQVver(fSv)  ( 

,  Smag )  1 

(2.56b) 

KMdispMTFVsys(fs  J  j 

CTF  (  fsh 

MdispMT  F Hsys  ( fSh ) 

(2.56c) 

CTF^U;-J 

MdispMT  F Vgys  (fSv ) 

(2.56d) 

The  Mdisp  is  the  display  glare  and  a  is  a  proportionality  constant  of  169.6  Hz  2. 
SCNtmp  is  solicited  from  the  user  and  considered  a  constant  throughout  calculations. 


5.  Detector  Cooling  in  NVThermlP 

NVThermlP  allows  for  the  selection  of  an  uncooled  detector.  Selecting  this 
option  requires  a  user  input  of  a  performance  measurement  for  the  uncooled  array 
in  terms  of  the  measured  detector  noise,  detector  frame  rate,  f-number,  and  optics 
transmission.  Then,  NVThermlP  will  calculate  a  peak  D*  and  integration  time  for 
the  uncooled  sensor. 
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6.  Targeting  Task  Performance  (TTP)  Metric  and  NV- 
ThermlP 

Unlike  FLIR92  and  previous  version  of  of  NVTherm,  NVThcrmlP  uses  the 
Targeting  Task  Performance  (TTP)  metric  to  predict  the  probability  of  successful 
task  performance.  Since  this  study  does  not  utilize  the  range  performance  predictions 
from  NVThermlP,  no  further  details  are  provided  here,  but  for  more  information,  see 
the  NVThermlP  User’s  Manual  (NVESD  2005). 
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III.  OPTICAL  TURBULENCE  MODEL 

CHOICE 

Optical  systems  operating  in  a  turbulent  atmosphere  experience  broadening  of 
the  point  spread  function,  which  is  best  described  using  a  specific  atmospheric  MTF. 
For  a  collimated  beam  passing  through  the  atmosphere,  turbulence  distorts  the  shape 
of  the  wavefront  and  causes  variations  of  intensity  along  the  wavefront  so  that  when 
the  beam  is  brought  back  into  focus  by  some  optical  system,  the  image  formed  has 
been  altered  by  the  atmospheric  affects.  Driggers  (1999)  notes  that  the  atmospheric 
MTF  may  be  expressed  as  two  separate  MTFs  created  by  two  different  effects  in  the 
atmosphere:  scattering  due  to  aerosols  and  blurring  due  to  turbulence.  In  this  study, 
the  MTF  describing  scattering  due  to  aerosols  is  neglected. 

This  study  uses  the  mathematical  model  of  the  atmospheric  turbulence  MTF 
developed  by  Goodman  (1985).  Various  mathematical  determinations  of  the  atmo¬ 
spheric  parameter  C%  are  considered  for  use  in  the  atmospheric  turbulence  MTF 
calculation,  in  order  to  most  appropriately  represent  C\  through  the  whole  atmo¬ 
sphere. 

A.  OPTICAL  TURBULENCE  AND  ATMOSPHERIC  C% 

According  to  the  IREO  Handbook  Vol.2  (1993),  turbulence  in  the  atmosphere 
creates  random  variations  in  the  atmosphere’s  index  of  refraction.  These  irregularities 
distort  the  wavefronts  that  pass  through  them  and  thereby  cause  image  distortion 
or  blur  in  imaging  systems.  Although  there  are  different  geometries  for  describing 
turbulent  systems,  the  applicable  one  in  dealing  with  airborne  and  ground-based 
sensors  is  spherical  wave  propagation.  In  this  case,  the  propagating  light  comes  from 
sources  that  are  in  or  near  the  turbulence,  as  in  the  imaging  of  objects  in  the  turbulent 
atmosphere.  The  sensors  receiving  the  light  are  also  in  the  turbulent  atmosphere. 
Although  this  study  focuses  on  the  airborne  sensor  situation  where  a  target  on  the 
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surface  is  imaged  by  an  over- flying  aircraft,  it  is  important  to  note  that  it  is  the 
turbulence  at  the  sensor  that  most  strongly  affects  the  image  quality. 

Operationally,  the  best  estimate  of  atmospheric  turbulence  would  be  deter¬ 
mined  from  an  analysis  or  forecast,  but  when  a  local  estimate  is  not  available,  the 
atmospheric  turbulence  may  be  predicted  using  models  that  estimate  the  atmospheric 
structure  parameter  C2,  as  described  below.  The  atmospheric  structure  parameter  C2 
describes  atmospheric  variations  in  the  index  of  refraction.  Essentially,  fluctuations 
in  the  index  of  refraction  along  the  path  between  the  target  and  the  sensor  causes 
some  IR  radiation  to  be  randomly  bent  from  the  path,  resulting  in  a  blurred  image 
at  the  sensor.  The  more  turbulent  the  atmosphere,  the  greater  the  blurring.  C2  is 
best  defined  by  Equation  3.57, 

C2(x)  =  [n{x)  —  n(x  +  r)]r~%  (3.57) 

where  n(x)  is  the  index  of  refraction  at  a  point  x  in  the  atmosphere  and  n(x+r )  is  the 
index  of  refraction  some  distance  r  away  from  x.  The  over-bar  indicates  an  average 
over  the  representative  part  of  the  atmosphere  by  either  averaging  n  over  time  at  one 
location  or  space  at  a  very  short  time  (Goroch  1980). 

Atmospheric  variations  in  the  index  of  refraction  are  described  by  C2.  The 
generic  form  of  how  C\  varies  in  the  atmosphere  is  given  by  Friehe  (1977)  as 

C2  =  (79  x  1(T6  — J  [C2  +  0.113Ctq  +  3.2  x  KT3^)  (3.58) 

where  P  is  the  pressure  in  mb,  T  is  the  temperature  in  K,  is  the  temperature 
structure  function  parameter,  Ctq  the  temperature-water  vapor  parameter  and  Cq 
the  water  vapor  parameter,  all  in  K2/m~i .  The  IREO  Handbook  Vol.2  (1993)  states 
that  the  dry-air,  or  Cj*,  term  dominates  in  most  applications  since  generally  the  C^q 
and  Cq  terms  comprise  no  more  than  2%  of  the  total  C2.  Since  Cf.  is  typically  the 
dominant  term  in  C2,  in  this  study  it  will  be  considered  the  only  contributor.  All 
three  structure  parameters  vary  depending  on  location  in  the  boundary  layer  (Fairall 
1982).  Although  it  will  not  be  addressed  in  this  study,  it  should  be  noted  that  in  the 
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microwave  part  of  the  spectrum,  C2  depends  much  more  strongly  on  the  humidity 
and  Ctq  and  Cq  are  large  contributors  to  the  overall  C\  (Fante  1980).  Also,  in 
rare  situations  in  the  IR  part  of  the  spectrum,  Ctq  nray  be  more  significant  ( IREO 
Handbook  Vol.2  1993). 

1.  Turbulence  Near  the  Ground 

There  are  two  primary  ways  that  the  ground  impacts  the  movement  of  air  and 
causes  turbulence.  First,  the  free  air  stream  flowing  along  the  ground  experiences 
friction  due  to  surface  roughnesses,  which  causes  wind  shear.  Second,  the  ground 
may  serve  as  a  source  or  a  sink  for  thermal  energy  of  the  air.  Given  sunny  conditions 
during  the  day,  the  ground  will  be  a  source  of  heat  because  the  sun  warms  the  ground 
to  a  temperature  higher  than  the  air  above  it.  This  leads  to  thermal  convection 
and  instability.  At  night,  the  ground  will  act  as  a  heat  sink  by  radiative  cooling, 
resulting  in  a  ground  temperature  colder  than  the  air  above  it.  These  conditions 
are  considered  stable.  When  the  air  and  the  ground  are  at  the  same  temperature, 
atmospheric  conditions  are  considered  neutral.  These  extreme  fluctuations  in  ground 
temperature  are  observed  over  land,  but  over  the  ocean  the  temperature  difference 
between  night  and  day  is  much  less.  The  variation  in  temperature  differences  was  the 
primary  motivation  for  comparing  FLIR92  and  NVThermlP  over  desert  and  marine 
locations. 

Fairall  et.  al.  (1982)  propose  that  in  the  surface  layer  of  the  atmosphere,  the 
structure  function  parameters  have  the  forms  of 

C2t  =  T?Z-t/(0  (3.59a) 

C2q  =  QlZ-2sAf(£)  (3.59b) 

Ctq  =  ttqT^Q^Z  zy/AffJi)  (3.59c) 

where  T*  and  Q*  are  temperature  and  humidity  scaling  parameters  in  K  and  g/m3 

respectively,  /(£)  is  a  dimensionless  function,  and  Ttq  is  the  temperature-humidity 
correction  parameter.  Fairall  et.  al.  (1982)  gives  an  estimate  value  of  0.8  for  Ttq  in 
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unstable  conditions  such  as  those  being  studied  here.  A  is  a  constant,  taken  to  be 
0.6.  Z  is  the  height  above  the  surface  in  m. 

The  dimensionless  function  /(£)  is  connected  to  the  Richardson  number,  so  it 
accounts  for  varying  conditions  depending  on  the  stability.  Wyngaard  et.  al.  (1971) 
determined  the  empirical  form, 


/(£) 

=  4.9(1  -  7.0£)“i 

£  <  o 

(3.60a) 

/(£) 

=  4.9 

£  =  o 

(3.60b) 

/(£) 

=  4.9(1  + 2. 75£) 

£  >  o 

(3.60c) 

where  £  <  0  corresponds  to  unstable  conditions  (as  under  consideration  here),  £  =  0 
is  neutral  and  £  >  0  is  stable.  Wyngaard  originally  split  /(£)  into  two  segments  of 
£  <  0  and  £  >  0,  but  the  neutral  case  has  been  added  here  to  show  when  /(£)  is  a 
constant.  £  itself  is  defined  as  the  height  scaled  by  the  Monin-Obukhov  length  scale, 


z  KgZ  (t*  +  0.61Ta^) 

L=  7X 


(3.61) 


L  is  the  Monin-Obukhov  length  scale  and  is  defined  as  the  height  over  the  ground 
where  the  mechanically  produced  turbulence  from  vertical  shear  balances  the  dissi¬ 
pative  effect  of  negative  buoyancy.  In  other  words,  the  Richardson  number  equals 
unity.  L  may  be  expressed  as 


L 


T  u 1 


Kg  (T*  +  0.61Ta^) 


(3.62) 


where  k  is  the  unitless  von  Karman  constant  that  Fairall  et.  al.  approximate  to  0.35, 
g  is  the  gravitational  acceleration  on  earth  taken  as  9.8  m/s2,  Q *  is  the  humidity 
scaling  parameter  in  g/m3,  p  is  the  density  of  air  in  kg/m3,  and  rt*  is  the  frictional 
velocity  in  m/s.  In  this  case,  the  density  of  air  is  taken  as  1.3  kg/m3.  Ta  is  the 
temperature  of  the  ambient  air  in  the  region  of  interest,  in  K. 

If  C2n  is  taken  to  primarily  depend  on  C/  as  was  assumed  above,  then  Equations 
3.59b  and  3.59c  may  be  neglected  and  their  terms  in  Equation  3.58  dropped.  Davidson 
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et.  al.  (1978)  indicate  that  above  the  ocean  and  within  10m  of  the  ocean  surface, 
Cj<  varies  as  but  above  that  10m,  Cl  varies  as  Z~ §.  This  is  better  seen  by 

solving  Equation  3.59a  for  the  conditions  given  in  Equations  3.60,  as  shown  below 
where  again  £  =  j- 


Cl 

=  A.9T2Z-i  [1  -  7£P 

£<o 

(3.63a) 

Cl 

=  4.9  T^Z~i 

£  =  o 

(3.63b) 

Cl 

=  A.9T2Z~%  (1  +  2.75£) 

£>o 

(3.63c) 

Substituting  Equation  3.61into  3.63  for  the  three  cases  of  unstable,  neutral,  and  stable 
C\  respectively,  yields 


cl 

3  * 

(3.64a) 

cl 

=  4.9  TlZ~% 

(3.64b) 

cl 

~  13.5  T°;Z^L~l. 

(3.64c) 

In  unstable  conditions,  |£|  3>  1  (Equation  3.64a),  so  (1  —  7£)  3  can  be  approximated 
_  2 

as  (— 7£)  3.  In  the  stable  case,  |£|  <C  1  (Equation  3.64c),  so  (1  +  2.75£)  is  approx¬ 
imately  2.75£.  Other  authors  indicate  that  for  |£|  <C  1,  Cl  may  even  be  taken  as 
independent  of  Z . 

The  above  discussion  applies  over  the  ocean,  but  Hall  (1977)  observed  that 
in  unstable  conditions  near  land  surface,  C'2  still  decreases  with  height  by  Z~  s.  The 

4 

Z~  3  height  dependence  frequently  continues  beyond  10m  over  land,  but  is  still  only 
valid  in  the  surface  layer. 

2.  Turbulence  Above  the  Surface  Layer 

In  the  so-called  mixed  layer  above  the  surface  layer,  Fairall  et.  al.  (1982)  offer 
definitions  for  C'l,  Cq,  and  Ctq  where  in  all  cases  C2n  falls  off  as  a  function  of 
They  are  very  similar  to  Equation  3.60,  so  in  this  study,  Equation  3.60  is  used  beyond 
the  lowest  surface  layers. 
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Above  the  mixed  layer,  one  possible  model  of  C2  behavior  is  the  Hufnagel 
model.  Goodman  (1985)  offers  a  possible  analytic  approximation  to  how  C2  varies  in 
the  vertical  and  references  the  work  by  Hufnagel  and  Stanley  (1964).  This  approxi¬ 
mation  is 

/  y  \  10 

C2n{Z)  =  2.7  x  10"16  3<u2>  ( —  J  +  e~™  (3.65) 

where  Z  is  the  location  along  the  flight  path  in  km  and  <  u2  >  is  the  mean  value 
of  the  squared  wind  speed.  C2  is  again  in  m~i .  This  approximation  reasonably  rep¬ 
resents  the  average  variation  of  C2  in  the  middle  to  upper  atmosphere  but  poorly 
represents  the  boundary  layer.  Variations  of  C2  in  the  boundary  layer  differ  exten¬ 
sively  depending  on  the  stability  of  the  boundary  layer,  as  described  above.  Fante 
(1980)  and  other  sources  indicate  that  the  Hufnagel  approximation  agrees  fairly  well 
with  observation  starting  around  5  km  above  the  surface. 

3.  Turbulence  Through  Whole  Atmosphere 

In  this  study,  the  atmospheric  structure  parameter  C2  is  taken  to  follow  the 
bulk  method  below  200  m  and  the  Hufnagel  formulation  above  the  5  km  point.  It  is 
worth  noting  that  the  Hufnagel  formulation  is  an  early  numerical  model  of  the  upper 
atmospheric  structure  parameter  and  is  generally  only  considered  valid  for  above  5 
km  in  the  atmosphere.  For  simplicity  in  MatLab  programming,  between  200  m  along 
the  pathlcngth  and  5  km,  an  average  value  of  1  x  10~16  was  chosen.  It  is  near 
the  surface  that  C2  varies  the  most  and  is  the  largest,  so  the  actual  choice  for  C2 
above  the  surface  layer  was  of  less  concern  than  that  in  the  surface  layer.  The  C2 
calculation  is  summarized  in  Equation  3.66, 


c2 

=  BulkMethod 

Z  <  200m 

(3.66a) 

r2 

=  lx  10 

200m  <  Z  <  5000m 

(3.66b) 

r2 

n 

=  Hufnagel  Method 

Z  >  5000m. 

(3.66c) 

In  order  to  test  the  validity  of  these  assumptions,  the  bulk  method  and  the  Hufnagel 
formulation  were  plotted  for  several  temperature  differences.  In  Figure  14,  the  surface 
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temperature  is  taken  to  be  28°C  and  the  ambient  air  temperature  25°C.  Other  param¬ 
eters  are  given  in  Table  IX,  where  for  this  test  case,  generic  parameters  are  chosen. 
The  two  methods  intersect  below  200  m,  but  at  200  m  the  C%  values  are  on  the  order 
of  1  x  10  16  m  3,  which  agrees  fairly  well  with  results  from  Fairall  et.  ah  (1982)  and 
other  sources.  Values  differed  by  less  than  an  order  of  magnitude  using  the  Hufnagel 
model  at  5  km.  The  choice  of  a  constant  in  the  middle  atmosphere  simplified  the 
MatLab  programming,  but  the  small  scale  of  and  the  relatively  small  difference 
between  the  two  models  around  the  500  m  to  5000  m  layer  suggested  that  choosing  a 
constant  Cfx  was  reasonably  representative  of  the  atmospheric  conditions.  Note  that 
for  simplicity  in  the  MatLab  programming,  the  bulk  method,  which  describes  C%  in 
the  surface  layer,  is  taken  to  be  valid  up  to  200  m.  This  does  not  mean  that  the 
surface  layer  itself  extends  up  to  200  m.  Rather,  the  bulk  method  is  assumed  valid 
up  to  200  nr,  so  the  bulk  method  may  be  assumed  valid  beyond  the  surface  layer. 

Table  IX.  Parameters  for  Test  Cases. 

Generic  parameters  chosen  to  test  the  Bulk  and  Hufnagel  formulations  of  C%- 


Parameter 

Value 

Units 

Rel  Humidity  Ambient  Air 

0.8 

fraction 

Rcl  Humidity  Sfc  Air 

1.0 

fraction 

Wind  Speed 

7 

m/s 

Sfc  Pressure 

1012 

hPa 

Thermal  Sfc  Roughness 

0.02 

m 

Momentum  Sfc  Roughness 

0.02 

m 

Total  Pathlcngth 

10,000 

m 

After  comparing  the  bulk  method  and  the  Hufnagel  method  for  a  sample 
test  case,  the  two  methods  were  compared  for  conditions  at  White  Sands  and  Point 
Conception  (data  give  in  Chapter  IV).  As  indicated  in  Figures  15  and  16,  the  assumed 
constant  value  of  1  x  10-16  m-I  appears  fairly  representative  between  the  two  models. 
Also,  the  exponential  behavior  of  the  bulk  method  near  the  surface,  in  both  cases, 
indicates  how  much  stronger  turbulence  is  near  the  ground. 
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Note  that  while  the  surface  bulk  model  and  Hufnagel  model  were  used  this 
study,  there  are  numerous  other  models  of  the  atmospheric  structure  parameter.  It 
is  also  important  to  note  that  this  study  does  not  consider  possible  surface  inversions 
and  neglects  the  stable  case. 


B.  MTF  FOR  OPTICAL  TURBULENCE 

The  Goodman  (1985)  formulation  for  atmospheric  optical  turbulence  was  cho¬ 
sen  for  this  study.  It  can  be  expressed  for  the  long-exposure  and  short-exposure  case, 
although  the  focus  here  is  on  the  long-exposure  case. 

The  long-exposure  and  short-exposure  atmospheric  turbulence  MTFs  differ 


Comparing  C2  Models  Sfc  T  =  28°C,  Amb  T  =  25°  C 


10 

C*  (rrf2/3) 


Figure  14.  Comparison  of  the  Bulk  and  Hufnagel  Methods  with  Sfc  Temp  28°C. 
At  28°C,  the  two  methods  intersect  at  158  nr.  At  200  m,  the  bulk  method  has  a  C* 
value  of  approximately  1  x  10~16  m“3.  The  bulk  method  is  used  below  200  m. 
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substantially.  In  the  long-exposure  case,  the  exposure  time  is  such  that  the  phase  and 
log-amplitude  vary  over  the  course  of  the  exposure.  Thus  during  a  long  exposure,  the 
image  recorded  will  be  spread  due  to  random  variations  in  the  tilt  of  the  wavefront. 
To  consider  the  short-exposure  case,  a  random  factor  associated  with  the  wavefront 
tilt  is  extracted  from  the  MTF  before  the  average  is  taken.  This  is  not  done  for  the 
long-exposure  case.  For  the  extremely  short  exposures,  there  will  be  no  impact  from 
the  wavefront  tilt  and  it  is  ignored  in  determining  the  MTF  (Fried  1966). 

The  derivation  for  the  Goodman  atmospheric  MTF  is  summarized  in  Ap¬ 
pendix  A.  Using  the  generic  atmospheric  input  parameters  given  in  Tabic  IX,  the 
Goodman  long  exposure  atmospheric  turbulence  MTF  is  orders  of  magnitude  smaller 
than  the  other  MTFs  included  in  FLIR92.  In  Figure  17,  Goodman  MTF  is  subtracted 
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Figure  15.  Comparison  of  the  Bulk  and  Hufnagel  Methods  at  White  Sands. 
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from  unity  to  better  illustrate  the  form  of  the  MTF.  The  used  in  the  Goodman 
calculations  may  be  calculated  using  Equation  3.66. 

From  the  above  discussions,  it  appears  that  C\  is  most  significant  near  the 
surface.  In  one  illustration  of  this,  the  Goodman  atmospheric  MTF  is  examined 
for  a  variety  of  slant  angles  starting  with  a  path  perpendicular  to  the  surface,  and 
ending  with  a  virtually  horizontal  path  through  the  surface  layer,  parallel  to  the 
surface.  Figure  18  illustrates  the  Goodman  MTF  at  various  slant  angles,  for  a  constant 
pathlength  of  100  m.  Plotted  is  unity  minus  the  Goodman  MTF,  as  in  Figure  17,  to 
best  observe  this  small  scale  MTF.  Although  Figure  18  makes  it  clear  that  even  close 
to  parallel  to  the  surface,  the  Goodman  MTF  is  small  scale,  still  the  resultant  MTF 
is  more  significant  at  smaller  slant  angles. 


Comparing  Models  Sfc  T  =  31C,  Amb  T  =  25°  C 


Figure  16.  Comparison  of  the  Bulk  and  Hufnagel  Methods  at  Pt  Conception. 
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Figure  17.  Goodman  Long  Exposure  Atmospheric  Turbulence  MTF. 

Note  the  small  scale  of  the  atmospheric  turbulence  MTF.  Plot  is  of  (1  -  Goodman 

MTF). 
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MTF 


Figure  18.  Comparison  of  the  Goodman  MTF  at  Different  Slant  Angles. 
Perpendicular  to  the  surface,  the  Goodman  MTF  is  least  significant  and  parallel  to 
the  surface,  most  significant.  In  all  cases,  the  Goodman  MTF  is  small  scale. 


IV. 


METHODS  AND  DATA  SETS 


A.  SELECTION  OF  DATA  LOCATIONS 

This  study  compares  MRTDs  for  marine  and  desert  environments,  with  and 
without  an  included  MTF  for  atmospheric  turbulence.  Since  this  study  is  only  an 
initial  comparison,  only  two  sites  were  selected  for  analysis.  Because  of  access  to  ac¬ 
curate  data,  the  White  Sands  Missile  Range  in  New  Mexico  was  selected  for  the  desert 
environment  and  Buoy  46063,  just  off  Point  Conception  near  Santa  Barbara,  CA  was 
selected  for  the  marine  environment.  Although  the  two  locations  differ  in  nearly  all 
atmospheric  conditions,  the  critical  difference  is  between  the  surface  temperature  and 
the  ambient  air  temperature.  As  explained  in  Section  Ill. A,  the  temperature  struc¬ 
ture  parameter  is  the  largest  contributor  to  variations  in  the  atmospheric  structure 
parameter. 

B.  SELECTION  OF  ATMOSPHERIC  DATA 

Atmospheric  data  for  inclusion  in  TAWS  and  for  calculation  of  parameters 
in  FL1R92  and  NVThermlP  were  obtained  from  several  sources.  For  the  land  lo¬ 
cations,  atmospheric  data  was  obtained  from  the  Weather  Underground  website, 
www.wunderground.com.  The  website  provided  all  the  needed  TAWS  inputs.  An 
extreme  example  of  the  surface  temperature  and  ambient  air  temperature  at  White 
Sands  was  provided  by  Rick  Sharkey  and  cohorts  at  the  White  Sands  facility.  These 
input  parameters  are  given  in  Table  X. 

For  the  atmospheric  data  over  marine  locations,  observations  from  ships  and 
buoys  were  used.  The  team  at  AFWA  provided  ship  data  and  upper  air  analysis  for 
the  Point  Conception  region.  The  input  parameters  are  given  in  Table  XI 

In  both  cases,  data  was  obtained  for  27  Aug  07,  at  2030  GMT. 
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Table  X.  Atmospheric  Parameters  around  White  Sands. 


Atmospheric  Parameter 

Value 

Units 

Ambient  Air  Temperature 

35 

°C 

Surface  Air  Temperature 

75 

°C 

Ambient  Relative  Humidity 

.33 

fraction 

Surface  Relative  Humidity 

1.0 

fraction 

Ground  Type 

sand 

- 

Wind  Speed 

3.1 

m/s 

Wind  Direction 

27° 

NNE 

Surface  Pressure 

1003 

hPa 

Visibility 

16 

km 

Lower  Level  Clouds 

2000 

m 

Lower  Level  Cloud  Type 

Scat 

- 

Mid  Level  Clouds 

none 

- 

Upper  Level  Clouds 

3700 

m 

Upper  Level  Cloud  Type 

Scat 

- 

C.  SELECTION  OF  TARGET  LOCATIONS 

In  each  case  studied,  the  target  location  entered  into  TAWS  corresponded  to 
a  coordinates  of  a  nearby  region.  These  coordinates  are  summarized  in  Table  XII. 


D.  PROCEDURES 

In  this  study,  the  MRTD  was  calculated  for  chosen  atmospheric  conditions 
using  both  FLIR92  and  NVThermIP. 

1.  Parameters  Used  in  Calculations 

The  parameters  used  for  the  FLIR92  and  NVThermIP  calculations  are  listed 
in  tables  in  Appendix  D.  Whenever  possible,  the  same  inputs  were  used  in  FLIR92 
and  NVThermIP. 

2.  FLIR92  Calculations 

The  FLIR92  calculations  were  two-fold.  First,  a  MatLab  version  of  the  cur¬ 
rent  operational  FLIR92,  supplied  by  Dr.  Andy  Goroch,  was  used.  Then  the  MatLab 
version  of  FLIR92  was  modified  to  include  atmospheric  turbulence.  In  order  to  deter- 
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Table  XI.  Atmospheric  Parameters  around  Point  Conception. 


Atmospheric  Parameter 

Value 

Units 

Ambient  Air  Temperature 

25 

°C 

Surface  Air  Temperature 

31 

°C 

Ambient  Relative  Humidity 

.55 

fraction 

Surface  Relative  Humidity 

1.0 

fraction 

Ground  Type 

water 

- 

Wind  Speed 

3.7 

m/s 

Wind  Direction 

290° 

WNW 

Surface  Pressure 

1016 

hPa 

Visibility 

16 

km 

Lower  Level  Clouds 

none 

- 

Mid  Level  Clouds 

none 

- 

LIpper  Level  Clouds 

none 

- 

Table  XII.  Target  Locations. 

The  target  coordinates  are  entered  into  TAWS. 


Location 

Coordinates 

N  White  Sands  Range 

33°47T0”N;106°27,47”W 

Buoy  46063,  Pt  Conception 

34o16’21”N;120o41’55”W 

mine  the  Monin- Obukhov  length  and  T*,  g*,  and  ?!*,  certain  atmospheric  conditions 
were  entered  into  the  MatLab  FLIR92  program,  as  outlined  in  Table  XIII  below. 
Most  actual  values  depended  on  the  atmospheric  conditions  and  are  listed  later  for 
the  specific  cases. 

Table  XIII.  Atmospheric  Condition  Variables  for  Input  into  Modified  FLIR92. 


Variable 

Units 

Air  Temp 

°C 

Sfc  Temp 

°C 

Rcl  Humidity  Air 

fraction 

Rcl  Humidity  Sfc 

fraction 

Wind  Speed 

m/s 

Sfc  Pressure 

hPa 

Thermal  Sfc  Roughness 

0.2  m 

Momentum  Sfc  Roughness 

0.2  m 
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3.  NVThermlP  Calculations 

No  modifications  to  NVThermlP  were  made,  but  the  output  was  modified  as 
discussed  in  Section  1.1. C. 

4.  TAWS  Calculations 

The  resultant  MRTD  curves  were  entered  into  TAWS  version  5.2  and  the 
maximum  range  detection  results  from  TAWS  were  compared.  The  target  used  in 
TAWS  was  the  155  mm  M795,  a  projectile  in  a  black  plastic  bag.  The  sensor  used  in 
TAWS  was  the  generic  1000  IR  sensor,  modified  with  the  MRTD  outputs  from  either 
FL1R92  or  NVThermlP.  Weather  inputs  are  given  in  Tables  X  and  XI.  The  test  cases 
were  run  at  a  height  of  100  m. 

For  the  MRTD  input,  a  table  of  spatial  frequency  and  corresponding  MRTD 
result  was  inserted  into  TAWS.  The  FL1R92  output  had  a  larger  frequency  interval 
than  NVThermlP.  To  account  for  this,  the  NVThermlP  MRTD  results  were  inter¬ 
polated  onto  the  FL1R92  results.  Thus,  the  same  number  of  frequency-MRTD  pairs 
were  entered  for  both  the  NVThermlP  and  FL1R92  cases. 
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V. 


RESULTS 


A.  ATMOSPHERIC  TURBULENCE  IN  FLIR92 

The  modified  FLIR92  including  an  MTF  for  atmospheric  turbulence  with  a 
varying  dependence  on  C%,  was  compared  to  the  original  FLIR92.  The  inputs  are 
given  in  Appendix  D  with  the  atmospheric  inputs  for  the  modified  FL1R92  given  in 
Section  1V.B  for  White  Sands,  NM.  The  resultant  MRTDs  and  system  MTFs  were 
compared  in  the  horizontal  and  vertical  directions,  as  shown  in  Tables  XIV  and  XV. 
M  indicates  the  modified  FLIR92  with  atmospheric  turbulence,  O  the  original.  fs  is 
in  mrad-1. 

Table  XIV.  MRTD  for  Modified,  Original  FLIR92  in  the  horizontal  direction. 
The  M  indicates  the  modified  FLIR92  while  the  0  indicates  the  original  FLIR92.  fs 
is  in  mrad-1  and  the  MRTD  is  in  K . 


fs 

MRTDm 

MTFsySM 

MRTD0 

MTFsyso 

0.9935 

0.0013 

0.8214 

0.0013 

0.8214 

2.9804 

0.0054 

0.4874 

0.0054 

0.4875 

4.9674 

0.0155 

0.2429 

0.0155 

0.2429 

6.9543 

0.0453 

0.1008 

0.0453 

0.1008 

8.9413 

0.1511 

0.0343 

0.151 

0.0343 

10.9283 

0.6072 

0.0094 

0.6069 

0.0094 

12.9152 

3.13 

0.0019 

3.1284 

0.002 

14.9022 

23.2345 

0.0003 

23.2193 

0.0003 

Figure  19  illustrates  how  the  difference  between  the  modified  and  original 
FL1R92  is  only  apparent  at  higher  frequencies.  Even  so,  the  scale  indicates  how  small 
is  the  impact  of  including  atmospheric  turbulence  in  the  MRTD  calculation. 

There  is  a  difference  in  resultant  MRTD  when  including  the  atmospheric  tur¬ 
bulence  MTF,  as  seen  at  the  higher  frequencies,  so  starting  at  10.9283  cycles/mrad 
in  the  horizontal  case.  The  difference  is  minimal,  though,  less  than  0.2%  of  the  hor¬ 
izontal  MRTD  and  vertical  MRTD.  The  difference  is  too  small  to  impact  the  range 
calculation  in  TAWS  since  the  input  for  MRTD  in  TAWS  is  only  of  accuracy  10-3  K. 
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Table  XV.  MRTD  for  Modified,  Original  FLIR92  in  the  vertical  direction. 

The  M  indicates  the  modified  FLIR92  while  the  O  indicates  the  original  FLIR92.  fs 
is  in  mrad-1  and  the  MRTD  is  in  K . 


fs 

MRTDm 

MTFsySM 

MRTD0 

MTFsyso 

0.9935 

0.0013 

0.8254 

0.0013 

0.8254 

2.9804 

0.0053 

0.5094 

0.0053 

0.5095 

4.9674 

0.014 

0.2751 

0.014 

0.2751 

6.9543 

0.0365 

0.1295 

0.0364 

0.1295 

8.9413 

0.1023 

0.0528 

0.1023 

0.0528 

10.9283 

0.3224 

0.0185 

0.3222 

0.0185 

12.9152 

1.1786 

0.0054 

1.178 

0.0054 

14.9022 

5.2264 

0.0013 

5.223 

0.0013 

16.8891 

30.9544 

0.0002 

30.9296 

0.0002 

Again,  Figure  20  indicates  the  differences  at  higher  frequencies  and  the  small  scale  of 
the  differences. 


Horiz:  Modified  FLIR92  -  Original  FLIR92  MRTD 


Figure  19.  Horiz  Difference  Between  MRTD  Curves  for  Modified,  Original  FLIR92. 
Differences  are  not  observed  until  higher  frequencies. 
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B.  COMPARING  FLIR92  AND  NVTHERMIP 
1.  Comparing  Resultant  MRTDs 

The  modified  FLIR92  MRTD,  including  an  MTF  for  atmospheric  turbulence 
with  a  varying  dependence  on  C%,  was  compared  to  NVThermIP.  As  explained  in 
Section  1.1. C,  the  CTF  results  from  NVThermIP  were  used  to  find  a  comparable 
MRTD.  From  both  FLIR92  and  NVThermIP,  the  MRTD  results  were  split  into  hori¬ 
zontal  and  vertical  MRTDs.  Comparisons  of  the  horizontal  MRTD  show  fairly  similar 
results,  as  seen  in  Figure  21.  The  downturn  on  the  FLIR92  MRTD  curve  is  due  to 
insufficient  eye  filtering.  NVThermIP  attempts  to  compensate  for  this  weakness,  and 
the  difference  is  apparent  in  the  upturning  curve  at  low  frequencies.  Overall,  the 
NVThermIP  MRTD  is  greater  because  NVThermIP  includes  an  improved  and  less 
pessimistic  calculation  to  incorporate  the  human  visual  response  to  the  display  (see 
Section  B.l.n). 

Comparing  the  vertical  MRTD,  as  in  Figure  22,  shows  that  once  again  at  low 


Vert:  Modified  FLIR92  -  Original  FLIR92  MRTD 


Figure  20.  Vert  Difference  Between  MRTD  Curves  for  Modified,  Original  FLIR92. 
Differences  are  not  observed  until  higher  frequencies. 
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spatial  frequencies,  FLIR92  turns  down  while  NVThermlP  turns  up.  At  high  spatial 
frequencies,  though,  the  two  curves  appear  to  agree  closely  in  terms  of  slope  and 
value.  Because  this  study  focuses  on  differences  between  the  two  MRTD  models  and 
the  horizontal  MRTD  exhibits  a  greater  difference,  the  horizontal  MRTD  is  used  as 
the  input  into  TAWS. 

2.  Comparing  Detection  Range 

The  MRTDs  were  calculated  with  FLIR92  and  NVThermlP  for  both  White 
Sands  and  Point  Conception. 

The  resultant  maximum  range  detection  from  TAWS  at  35%,  65%,  and  95% 
probabilities,  are  compared  for  both  the  White  Sands  site  and  the  Point  Conception 
site.  NVThermlP  uses  a  constant  atmospheric  structure  parameter.  C\  was  set 


Figure  21.  Horizontal  MRTD  Curves  for  FLIR92,  NVThermlP. 

The  horizontal  MRTD  curves  have  different  behaviors  at  the  low  frequencies  and 

high  frequencies. 
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to  ICC14  m~i  for  both  cases.  Other  inputs  to  TAWS  are  listed  in  Section  IV. B. 
The  White  Sands  case  shows  NVThcrmlP  maximum  detection  range  going  from  0% 
to  1.5%  greater  than  the  FLIR92  maximum  detection  range,  for  all  probabilities. 
In  Figure  23,  the  maximum  detection  ranges  at  all  probabilities  for  MRTD  from 
NVThermlP  (top)  and  FLIR92  (bottom)  are  shown. 

Figure  24  better  illustrates  the  differences  between  the  maximum  detection 
ranges  calculated  using  the  MRTDs  from  NVThermlP  and  FLIR92.  The  largest 
percent  difference  is  at  35%  probability  of  detection. 

The  Point  Conception  case  shows  NVThermlP  maximum  detection  range  go¬ 
ing  from  0%  to  2%  greater  than  the  FLIR92  maximum  detection  range,  for  all  prob¬ 
abilities.  Figure  25  illustrates  the  maximum  detection  ranges  for  35%,  65%,  and  95% 
probabilities  of  detection  based  on  the  MRTDs  from  both  NVThermlP  (top)  and 
FLIR92  (bottom). 


Figure  22.  Vertical  MRTD  Curves  for  FLIR92,  NVThermlP. 
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To  better  examine  how  the  maximum  detection  ranges  based  on  MRTDs  from 
FLIR92  and  NVThermlP  compare,  Figure  26  shows  FLIR29  and  NVThermlP  max¬ 
imum  detection  range  results  for  each  probability.  The  largest  percent  difference  is 
again  at  35%  probability  of  detection. 

3.  Comparing  Marine/Desert  Conditions 

After  comparing  maximum  detection  ranges  based  on  MRTD  calculations  from 
FLIR92  and  NVThermlP,  differences  between  maximum  detection  ranges  at  White 
Sands  and  Point  Conception  were  compared.  As  shown  in  Figure  27,  at  35%  prob¬ 
ability  of  detection,  there  is  a  greater  difference  between  maximum  detection  range 
from  the  FLIR92  MRTD  and  the  NVThermlP  MRTD  for  Pt  Conception.  At  95% 
probability  of  detection,  as  shown  in  Figure  27,  there  is  much  less  difference  between 
maximum  detection  ranges  using  FLIR92  and  NVThermlP  MRTDs,  but  there  is  still 
an  observeablc  difference. 

The  difference  between  NVThermlP  and  FLIR92  at  each  location  is  better 
seen  in  the  bar  plots  in  Figures  28  and  29  showing  the  viewing  angle  versus  the 
maximum  detection  range  at  35%  and  95%  probabilities  of  detection  respectively. 
The  35%  case  most  clearly  demonstrates  that  at  White  Sands,  over  the  hot  desert, 
there  is  the  greatest  difference  in  the  maximum  detection  range  between  FLIR92 
and  NVThermlP,  with  NVThermlP  having  a  greater  range.  At  95%  probability  of 
detection,  the  difference  between  the  desert  and  the  marine  locations  is  much  smaller. 


Figure  23.  Comparing  Maximum  Detection  Range  at  White  Sands  using  FLIR92, 
NVThermIP. 

Polar  Plot  of  viewing  angle  (deg)  versus  maximum  detection  range  (km)  for  35%, 

65%,  and  95%  probabilities  of  detection,  using  the  MRTD  from  NVThermIP  (top) 

and  FL1R92  (bottom). 
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Max  Det  Rng  35%,  65%,  95%  Prob  White  Sands 


WS  Comparing  Max  Detection  35%  Prob 


Comparing  Max  Detection  65%  Prob 


180 


180  4 


Comparing  Max  Detection  95%  Prob 
180  4 


- NVThermlP 

—  FLIR92 


Figure  24.  Maximum  Detection  Range  for  White  Sands  at  Different  Probabilities. 
Polar  Plot  of  viewing  angle  (deg)  and  maximum  detection  range  (km)  for  both 


FLIR92  and  NVThermlP. 
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NVThermlP  Max  Det  Rng  35%,  65%,  95%  Prob,  Pt  Conception 
180  , 


- 65% 

■  ■  ■  ■  1  95%  1 


FLIR92  Max  Det  Rng  35%,  65%,  95%  Prob,  Pt  Conception 
180  5 


- 65% 

- 95% 


Figure  25.  Comparing  Maximum  Detection  Range  at  Pt  Conception  using  FLIR92, 
NVThermlP. 

Polar  Plot  showing  view  angle  (deg)  versus  the  maximum  detection  range  (km)  for 
35%,  65%,  and  95%  probabilities  of  detection  using  MRTD  from  NVThermlP  (top) 

and  FLIR92  (bottom). 
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Max  Det  Rng  35%,  65%,  95%  Prob,  Pt  Conception 


PC  Comparing  Max  Detection  35%  Prob 
180  ,, 


Comparing  Max  Detection  65%  Prob 
180  , 


Comparing  Max  Detection  95%  Prob 
180  4 


Figure  26.  Maximum  Detection  Range  for  Pt  Conception  at  Different  Probabilities. 
Polar  Plot  showing  view  angle  (deg)  versus  the  maximum  detection  range  (km)  for 


both  FLIR92  and  NVThermIP. 
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Figure  27.  Comparing  Maximum  Detection  Ranges  at  White  Sands  and  Pt  Concep¬ 
tion  at  35%  Prob  and  95%  Prob. 

Polar  Plot  showing  view  angle  (deg)  versus  the  maximum  detection  range  (km)  for 
both  FLIR92  and  NVThermlP  at  both  White  Sands  (WS)  and  Pt  Conception  (PC) 
at  35%  and  95%  probabilities  of  detection. 
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White  Sands:  NVThermlP  and  FLIR92  Max  Det  Rng  35% 


Pt  Conception:  NVThermlP  and  FLIR92  Max  Det  Rng  35% 


Figure  28.  Comparing  Maximum  Detection  Ranges  at  White  Sands  and  Pt  Concep¬ 
tion  at  35%  Prob. 

Bar  plot  showing  view  angle  (deg)  versus  the  maximum  detection  range  (km)  for 
both  FLIR92  and  NVThermlP  at  White  Sands  and  Pt  Conception  at  35% 

probability  of  detection. 
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White  Sands:  NVThermlP  and  FLIR92  Max  Det  Rng  95% 


Viewing  Angle  (deg) 


Pt  Conception:  NVThermlP  and  FLIR92  Max  Det  Rng  95% 


Viewing  Angle  (deg) 


Figure  29.  Comparing  Maximum  Detection  Ranges  at  White  Sands  and  Point  Con¬ 
ception  at  95%  Prob. 

Bar  plot  showing  view  angle  (deg)  versus  the  maximum  detection  range  (km)  for 
both  FLIR92  and  NVThermlP  at  White  Sands  and  Pt  Conception  at  95% 

probability  of  detection. 
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VI.  CONCLUSIONS  AND 
RECOMMENDATIONS 

This  study  considered  two  cases  to  test  possible  improvements  to  the  MRTD 
input  into  TAWS.  First,  FLIR92  was  modified  to  include  atmospheric  turbulence  and 
the  resultant  MRTDs  from  the  original  and  modified  FLIR92  were  compared.  Initial 
results  indicate  that  while  there  is  a  difference  in  including  atmospheric  turbulence 
in  the  MRTD  calculation,  the  difference  in  most  cases  is  too  small  to  be  significant. 
Second,  FLIR92  and  NVThermlP  MRTD  results  were  entered  into  TAWS  and  the 
resultant  maximum  detection  ranges  were  compared.  Again,  the  differences  in  the 
representative  test  cases  were  minimal.  In  comparing  the  FLIR92  and  NVThermlP 
results  in  TAWS,  the  percentage  improvement  in  using  NVThermlP  was  slightly 
greater  for  the  marine  test  case  than  for  the  desert  test  case. 

A.  CONCLUSIONS  FLIR92  WITH/ WITHOUT  ATMO¬ 
SPHERIC  TURBULENCE 

The  results  using  a  test  case  as  described  in  Chapter  IV  indicates  very  little 
difference  in  the  MRTD  when  an  atmospheric  turbulence  MTF  is  included  in  the  total 
system  MTF.  The  atmospheric  turbulence  MTF  is  very  close  to  unity,  on  the  order 
of  1  —  ICC7,  so  the  total  system  MTF  is  simply  being  multiplied  by  a  factor  of  unity. 
In  the  Goodman  long  exposure  MTF  given  in  Equation  A. 103,  the  general  form  is  an 
exponential  that  depends  on  fs,  A,  and  C%(z)  where  z  is  the  height  in  the  atmosphere. 
Three  possible  cases  exist  where  the  atmospheric  MTF  may  become  significant.  First, 
if  the  spatial  frequency  is  very  large,  corresponding  to  imaging  very  small  targets,  the 
MTF  may  be  significant.  Second,  if  the  diffraction  wavelength  is  very  small,  well 
below  the  bounds  of  an  IR  study,  the  MTF  may  be  significant.  Finally,  if  is  large, 
at  a  minimum  on  the  order  of  10~10  m-i ,  then  the  atmospheric  MTF  may  become 
more  important. 
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Near  the  surface  where  the  bulk  model  is  applicable,  C2  depends  on  the  am¬ 
bient  air  and  surface  temperatures,  wind  speed,  surface  pressure,  relative  humidity, 

4 

thermal  surface  roughness  and  momentum  surface  roughness  as  well  as  z~a.  Figure 
30  shows  that  for  the  White  Sands  case  near  the  surface,  C%  is  slightly  greater  than 
10~10  m-§,  indicating  that  only  very  close  to  the  surface  may  the  atmospheric  MTF 
be  important.  Too  close  to  the  surface,  though,  at  less  than  0.2  m,  the  bulk  model 
for  C2n  including  approximations  for  T*,  w*,  and  g*  may  no  longer  be  valid  due  to  the 
extreme  heat  fluxes  and  bumps  along  the  surface  itself  that  are  likely  greater  than 
the  sensing  height  and  will  affect  the  scaling  parameters. 


Bulk  Model  Near  Sfc  T  ,  =  35°C,  T.  =  75°C 

amb  sfc 


Figure  30.  Near  Surface  Behavior  of  the  Bulk  Model. 
The  height  above  the  surface  shown  here  is  0  to  1  /mi. 


B.  CONCLUSIONS  FLIR92  VERSUS  NVTHERMIP 

Comparing  the  maximum  detection  range  using  the  horizontal  MRTD  from 
FLIR92  and  NVThermlP,  there  is  not  a  significant  difference  between  the  results 
in  the  marine  environment  and  the  results  in  a  desert  environment.  In  both  cases, 
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the  MRTD  from  NVThermlP  yielded  a  greater  maximum  detection  range  by  up  to 
2%,  but  the  results  did  not  even  consistently  improve  by  2%.  Figures  21  and  22 
show  substantial  difference  between  FLIR92  and  NVThermlP  MRTD  results,  but 
apparently  these  differences  do  not  significantly  impact  TAWS  calculations. 

C.  RECOMMENDATIONS 

Since  neither  the  modification  of  FLIR92  nor  the  use  of  NVThermlP  instead 
of  FLIR92  resulted  in  significant  differences  in  TAWS,  it  seems  best  at  the  moment 
to  continue  using  FLIR92  to  calculate  the  MRTD  for  TAWS.  However,  should  TAWS 
increase  its  sensitivity  to  the  MRTD,  the  differences  described  in  this  study  may 
become  significant. 

Including  an  atmospheric  turbulence  MTF  in  FLIR92  simply  involves  multi¬ 
plying  the  total  system  MTF  by  the  atmospheric  MTF.  Such  a  simple  adjustment 
may  prove  useful  for  cases  where  C%  is  large.  The  same  adjustment  could  be  made 
in  NVThermlP,  where  the  current  atmospheric  MTF  assumes  a  constant  C\.  In  NV¬ 
ThermlP,  the  maximum  C\  is  10-10  m“§  but  using  such  a  large  up  through  the 
atmosphere  is  highly  unrealistic. 

If  the  use  of  FLIR92  for  MRTD  calculations  for  TAWS  continues,  other  mod¬ 
ifications  will  need  to  be  made.  FLIR92  assumes  a  CRT  display,  but  in  many  ap¬ 
plications  now  the  images  are  displayed  on  a  flat  panel.  It  is  possible  to  enter  a 
user-defined  display  MTF  describing  the  flat  panel,  but  updating  the  program  would 
be  a  better  solution. 
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APPENDIX  A.  GOODMAN  DERIVATION  FOR 

LONG  EXPOSURE 

The  following  derivation  is  primarily  summarized  from  Goodman  (1985),  with 
substantiations  from  other  sources. 

A  generally  accepted  form  of  the  refractive  index  of  air  at  optical  frequencies 
is  given  by  Goodman  as 

p 

n  =  1  +  77.6(1  +  7.52  x  1(T3A“2)  —  x  10"6  (A.67) 

where  A  is  the  wavelength  of  light  in  (m i,  P  is  the  atmospheric  pressure  in  mb,  and 
T  is  the  atmospheric  temperature  in  K.  To  simplify,  consider  that  pressure  variations 
are  small  compared  to  the  temperature  fluctuations  at  optical  frequencies  and  so  may 
be  neglected.  Then,  choosing  a  representative  A  =  0.5/jrn,  it  is  apparent  that  for  an 
incremental  change  dT  of  temperature, 

7Q  p 

dn  =  — —  x  10-6dT.  (A. 68) 

Inhomogeneities  in  the  refractive  index,  otherwise  described  as  turbulent  ed¬ 
dies,  can  be  considered  individual  packets  of  air  with  individual  characteristic  refrac¬ 
tive  indices.  Then,  the  power  spectral  density  of  homogeneous  turbulence,  <!>„(/?), 
where  k  is  the  vector  wavenumber  k  =  (kx,  fty,  Hz),  yields  a  measure  of  the  relative 
number  of  eddies  of  size  Lv  =  —  and  LY  =  —  and  Lz  =  — .  For  isotropic  turbu- 
lence,  $re(/c)  is  only  a  function  of  the  non- vector  wavenumber  k  that  corresponds  to 
an  eddy  of  size  L  where  L  —  — . 

In  the  classic  theory  of  turbulence  outlined  by  Kolmogorov  (1961)  and  ex¬ 
panded  by  Goodman,  the  power  spectral  density  consistss  of  three  different  regions 
depending  on  the  size  of  k.  Kolmogorov  hypothesizes  that  the  kinetic  energy  associ¬ 
ated  with  the  larger  eddies  will  be  redistributed  to  successively  smaller  and  smaller 
eddies  without  loss  of  energy.  For  the  smallest  k,  corresponding  to  large  scales  since 
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k  =  1/A,  turbulence  depends  on  large-scale  topographic  and  meteorological  condi¬ 
tions.  This  makes  it  unlikely  that  the  turbulence  would  be  isotropic  or  homogeneous, 
so  no  mathematical  form  of  turbulence  is  predicted  in  this  region.  Once  k  grows 
beyond  some  critical  wave  number  kq,  $n(/c)  can  be  determined  by  physical  laws 
governing  how  the  large  turbulent  eddies  break  up  into  smaller  ones.  In  this  re¬ 
gion  beyond  k0,  the  so-called  inertial  subrange  of  the  spectrum,  the  power  density 
spectrum  may  be  written  as 

<MK)  =  (A-69) 

where  again  C2  is  the  atmospheric  structure  constant.  Goodman  further  defines 
C2n  as  the  structure  constant  of  the  refractive  index  fluctuations  and  a  measure  of 
the  strength  of  the  fluctuations.  In  the  third  region,  beyond  another  critical  wave 
number  Krn ,  turbulent  eddies  are  very  small  scale  and  viscous  forces  may  cause  them 
to  dissipate  their  energy.  This  inner  scale  of  the  turbulence  results  in  rapid  decay  of 
<hn  where  k  >  nm.  The  decay  can  be  modeled  as 


‘M'O  =  3  exP  (~^)-  (A-7°) 

Tatarski  (1961)  points  out  that  this  model  is  only  valid  when  for  the  inner  scale  of 
the  turbulence,  /0  —  ,  Km.  is  chosen  equal  to  and  k  >  n0.  The  above  expression 

approaches  infinity  at  the  origin,  so  Goodman  uses  the  von  Karman  spectrum  instead. 
This  form  also  accounts  for  the  finite  amount  of  air  in  the  atmosphere  and  does  not 
become  arbitrarily  large  as  k  goes  to  0.  Approximated,  the  von  Karman  spectrum 
yields  a  power  spectrum  density  of 


30  (hi2  +  Kq 


(A. 71) 


and  Goodman  emphasizes  that  this  mathematical  form  is  only  an  artificial  means  of 
avoiding  the  poles.  The  von  Karman  spectrum  depends  on  the  critical  wave  numbers 
in  both  the  medium  and  small  size  turbulent  eddies.  Figure  31  illustrates  the  form 
of  the  von  Karman  spectrum  compared  to  the  Kolmogorov  spectrum. 
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Kolmogorov  vs  von  Karman  Spectrum 


Figure  31.  The  Kolmogorov  and  von  Karman  Spectra. 

The  Kolmogorov  spectrum  approaches  infinity  as  k  goes  to  zero,  but  the  von 
Karman  spectrum  avoids  the  division  by  zero  and  has  a  finite  value  at  k  =  0. 


The  structure  function  of  the  refractive  index  fluctuations  in  the  case  of  at¬ 
mospheric  turbulence  in  imaging  systems  is  defined  as 

Dn(fi,  r2)  =<  K(rl)  -  ni(r^)]2  >  (A.72) 

where  n\  is  the  random  refractive  index  and  f  =  f[  —  =  (Ax,  Ay,  A z).  Assuming 

homogeneous  turbulence, 

Dn(r)  —  2[Tn(0)  —  rn(r)]  (A.73) 

where  the  autocorrelation  function  Tn  may  be  expressed  as 

OO 

rn('^)  =  JJJ  ^n(«)  exp  (-jK  ■  r)d3K,  (A. 74) 

—  OO 

and  it  follows  that 

OO 

Tn(°)  =  IIJ  $n(Z)d3K.  (A. 75) 
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Combining,  then, 


OO 


Dn{r)  =  2  JJJ  [1  -  exp  (-jn  ■  r)]  <$>n(n)d3n.  (A.76) 

—  OO 

Assuming  a  real  solution  so  the  exponential  simplifies  to  cos(n-r),  and  when  Tn(f)  is 
spherically  symmetric  and  ri\  is  statistically  isotropic,  the  Fourier  transforms  of  the 
power  density  function  and  the  autocorrelation  function  may  be  expressed  as  single 
integrals.  In  this  case,  for  n  =  n2  +  n2  +  n2  and  r  =  \J (5x)2  +  (Sy)2  +  ( Sz )2, 

.  00 

47T  [' 

r n(r)  =  —  /  &n(K)Ksm(Kr)dK.  (A. 77) 

r  J 
0 

Then  the  structure  function  Dn  simplifies: 

OO 

Dn(r)  =  2rn(0)  —  8n  [  $„(«;)  — sin (nr)dn.  (A. 78) 

J  r 

0 

To  determining  Tn(0),  the  limit  as  r  — >  0  is  taken, 

00  •  /  \ 

lim  Tn(r)  =  lim47r  f  $n(n)K2'- - - — dn.  (A. 79) 

r^O  r^O  J  nr 

0 

Then,  since  the  integral  does  not  depend  on  r,  the  limit  may  be  moved  inside  the 


integral  to  find 


lim  Tn(r)  =  in  [  $n(«)«:2lim— — - — dn. 
r—>  0  J  r—>0  nr 


Taking  the  limit,  <Fr)(r)  as  r  approaches  0  is  found  to  be 

OO 

lim  Tn(r)  =  in  &n(n)n2dK 
r—>0  J 

0 

which  yields  a  Dn(r )  of 


Dn (r )  =  8nJ  ^n(n)n2 
0 


1  - 


sm  nr 


nr 


dn. 


(A. 80) 


(A.8i) 


(A. 82) 


Assuming  that  the  power  density  has  the  form  of  Equation  A. 69,  the  structure  func¬ 


tion  may  be  simplified  to 


Dn(r) 


8n 

30 


C; 


OO 

2  f  -5 
1  n  3 


0 


sin  (nr) 
nr 


dn. 


(A. 83) 
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Equation  A. 83  is  integrated  using  the  integral  identity 


x 


sin  (ax  J 
ax 


dx  = 


T(z/)  sin( 


■KV  N 
2  > 


ll+v 


where  —3  <  v  <  —1.  It  is  valid  to  use  Equation  A. 84  because  in  this  case  v 


D  (r)  =  — 
”  '  30 


5 


r  I  --  I  shi  (  J 


57 T\ 


C‘ir  3 


DJr )  =  lx  Clr  3  =  Clr* 


(A. 84) 


=  80 

(A. 85a) 
(A. 85b) 


where  the  coefficient  of  unity  in  Equation  A. 85b  is  due  to  the  chosen  definition  of  C^. 

Goodman’s  derivation  continues  as  follows  by  determining  the  optical  transfer 
function  (OTF).  The  derivation  will  be  summarized  here,  but  the  specific  details  are 
available  in  Goodman’s  text  (1985).  The  MTF  for  atmospheric  turbulence  is  simply 
the  modulus  of  the  OTF,  and  can  easily  be  calculated.  In  this  case,  the  OTF  is  a 
real  value,  so  the  modulus  of  the  OTF  is  the  same  as  the  OTF,  so  the  MTF  equals 
the  OTF.  This  derivation  focuses  on  long-exposure  systems  where  the  exposure  time 
is  much  greater  than  one  one-hundredth  of  a  second.  The  short-exposure  case  can 
also  be  significant  in  applications  but  is  neglected  in  this  study.  Also,  the  long- 
exposure  OTF  is  assumed  to  be  identical  to  the  ensemble-averaged  OTF.  Using  this 
assumption,  Goodman  derives  a  total  average  OTF  of  the  form 


OTF(v)  =  OTF0(is)OTFl(v ) 


(A. 86) 


where  v  =  \J ufj  +  Fy  and  ujj  and  vy  are  spatial  frequency  variables.  To  simplify  the 
derivation,  a  circularly  symmetric  OTF  is  assumed,  and  the  long-exposure  OTF  of 
the  atmosphere,  OTFl ,  has  the  form, 


OTFl  =  exp 


(A. 87) 


where  again  D(Xfu)  is  the  wave  structure  function,  A  =  |  where  c  is  the  speed  of 
light,  and  V  is  the  center  frequency  of  the  waveband  under  consideration. 
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In  the  near-field  calculation  of  the  wave  structure  function,  a  detailed  expres¬ 
sion  for  the  wave  structure  function  is  necessary  to  specifically  define  the  form  of  the 
long-exposure  atmospheric  OTF.  Three  assumptions  are  made. 

First,  the  target  is  assumed  to  be  very  far  away  from  the  detector  lens.  In 
this  way,  the  angular  extent  of  the  target  is  small  enough  that  over  the  long  time 
average,  the  entire  target  is  affected  uniformly  by  the  atmosphere.  This  assumption  is 
considered  the  ’isoplanatic’  assumption,  in  other  words  assuming  that  the  point  spread 
function  is  uniform  over  the  entire  FOV.  In  the  long-time-average,  this  assumption  is 
probably  a  close  approximation  to  truth  and  thus  is  not  very  restrictive. 

Second,  the  atmospheric  turbulence  is  assumed  to  be  homogeneous  and  isotropic 
over  a  finite  distance  z  in  front  of  the  detector  imaging  lens.  This  assumption  is  in¬ 
valid  when  considering  vertical  viewing,  which  is  a  problem  for  situations  such  as 
an  aircraft  imaging  a  ground-based  target.  This  assumption  will  be  removed  later 
though,  as  C„  varies  over  the  pathlength. 

Third,  the  imaging  system  is  assumed  to  be  so  far  within  the  near  field  of 
the  relevant  turbulent  eddies  that  each  ray  incident  on  the  inhomogeneous  medium 
may  be  approximated  to  merely  be  delayed  by  that  medium  with  no  bending  of  the 
rays.  In  other  words,  the  turbulence  is  considered  weak  enough  that  no  significant 

i 2 

amplitude  scintillation  effects  occur.  This  assumption  is  only  valid  when  z  <C  -f. 
Clearly,  in  actual  practice  this  assumption  is  also  unrealistic,  but  will  be  removed 
later. 

Without  going  into  the  details  of  the  derivation,  these  assumptions  are  used 
to  find  that  the  wave  structure  function  equals  the  phase  structure  function  and 
the  atmospheric  OTF  may  be  equivalently  expressed  in  Equation  A. 87  and  below  in 
Equation  A. 88, 

OTFl  =  exp  [— i/A^A/V)]  (A.88) 
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where  of  course  Ds(\fv)  is  the  phase  structure  function.  Using  the  covariance  func¬ 
tions  Cn,  not  to  be  confused  with  the  atmospheric  structure  constant  C%, 

Z  OO 

Ds[r )  =  k2  j  dz'  j  2 Cco  (z'  -  z")  -  2Cn\j (z'  -  z"f  +  r2  dz"  (A.89) 

0  0 

where  z  again  is  the  distance  from  the  target  to  the  imaging  lens.  The  difference  of  the 
covariance  functions  may  further  be  expressed  as  a  difference  of  structure  functions, 
Am 

£  2:  _ 

Ds{r )  =  k2  J  dz'  j  Dn\J (zf  —  z")2  +  r2  -  Dn  (, z '  -  z")  dz" .  (A.90) 

0  0 

Using  the  simplification  that  the  integrand  is  an  even  function  of  z'  —  z" ,  the  double 
integral  may  be  reduced  to  a  single  integral  such  that  with  Az  —  z!  —  z!' , 

Z 

Ds(r)  =  2k2J  (z  —  Az)  Dn\J (A z)2  +  r2Dn  (A z)  —  d(Az).  (A. 91) 

0 

Next,  the  specific  form  of  the  structure  function  that  was  found  above,  Equa¬ 
tion  A. 85b,  is  substituted  into  Equation  A. 91  to  yield  an  expression  for  the  phase 
structure  function, 

Z 

Ds(r)  =  2 k2Cl  J{z-  Az)  [([Az]2  +  r2)^  -  [Az}%]d(Az).  (A.92) 

0 

Because  Equation  A. 85b  is  only  valid  for  l0  <  r  <  L0,  Equation  A.92  is  only  valid  for 
A^  <  Lq.  This  makes  it  appear  that  Equation  A.92  is  only  valid  when  the  pathlength 
z  is  less  than  the  outer  scale  Lo,  but  in  the  typical  cases  where  the  path  separation 
r  is  less  than  A,  the  integrand  vanishes  for  large  A2  and  it  does  not  matter  what 
was  the  exact  form  of  the  structure  function  for  r  >  A-  For  Az  >  A  and  r  <C  A, 
the  integrand  approaches  zero.  Thus,  As(r )  can  be  used  for  all  pathlcngths  without 
violating  any  already  stated  assumptions. 

Since  Az  is  much  less  than  z,  Equation  A.92  may  be  rewritten,  neglecting  the 
term  of  Az  multiplied  by  the  rest  of  the  integrand, 

OO 

Ds{r)  ^  2k2C2z  J  [([Az]2  +  r2)^  -  [Az]§]d(Az)  (A.93) 

0 
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where  the  upper  integration  bound  has  been  extended  to  oo  since  from  z  to  oo  barely 
contributes  to  the  integrand.  This  integral  may  be  evaluated  numerically  using  a  pro¬ 
gram  such  as  Mathematica.  First,  a  change  of  variables  makes  the  integral  expression 
clearer:  Az  =  ru  and  dAz  =  rdu, 

OO 

Ds(r)  =  Ds{r)  =  k2C^z  J  ( u 2  +  l)3  —  du.  (A. 94) 

—  OO 

Mathematica  gives  the  result  that 

Ds(r )  =  2.9143^(7^1.  (A.95) 


Substituting  Equation  A.95  into  the  expression  for  the  long-exposure  atmospheric 
OTF,  Equation  A. 88,  the  final  expression  for  the  long  exposure  atmospheric  OTF  is, 

OTFl(u)  =  exp  —  ^  x  2.9lP(7^z  (A/v) 3  .  (A. 96) 

To  obtain  a  form  of  OTFl  that  is  independent  of  the  optical  system,  the  OTF 
is  expressed  in  terms  of  the  frequency  fs  measured  in  cycles/rad  rather  than  v  in 
cycles/m.  In  that  case,  fs  =  fv  and 

OTFlUs)  =  exp  [-1  x  2.9lPC^  (A/,)1]  .  (A.97) 

To  further  simplify,  k  =  so  Equation  A.97  becomes 

^2^2^  ft 

OTFL(fs)  =  exp  -5.82  -V Js 

X  3 

The  long  exposure  OTF  is  identical  to  the  long  exposure  MTF  since  the  MTF  is 
merely  the  modulus  of  the  OTF  and  in  this  case  the  modulus  is  equal  to  the  OTF 
itself.  Also,  though  the  above  result  is  limited  to  the  three  assumptions  stated  above, 
it  may  be  generalized. 

First,  the  assumption  that  refractive  index  perturbations  only  delayed  rays 
rather  than  bending  or  diffracting  them  may  be  eliminated  without  any  change  to 
the  long  exposure  OTF  theory.  Goodman  proceeds  through  an  analysis  considering 
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the  case  of  refractive  index  perturbations  occurring  within  a  finite  region  along  the 
propagation  path  and  then  integrates  along  the  entire  propagation  path.  This  gener¬ 
alization  leads  to  exactly  the  same  results  as  for  the  previous  analysis,  indicating  that 
it  is  possible  to  consider  cases  where  the  rays  diffract  and  bend  as  they  pass  through 
varying  refractive  indices. 

Second,  since  variations  in  Cl  occur  frequently  in  practice,  the  assumption  that 
the  turbulence  is  uniform  is  removed.  The  analysis  in  Goodman  examines  smooth 
and  slow  variations  to  the  atmospheric  structure  constant  C\  over  the  path  length  in 
both  the  vertical  and  horizontal  directions.  Goodman  points  out  that  his  analysis  is 
an  approximation,  but  has  been  verified  with  more  rigorous  methods.  Since  Good¬ 
man  assumes  a  slow  and  smooth  change  in  C%,  his  method  of  determining  the  OTF 
considers  a  ”  quasihomogeneous”  model  for  the  structure  function  of  the  refractive 
index  fluctuations, 

Dn(n,  r2)  =  C2n  \ri  ~  r2\%  (A.99) 

where  /0  <  r)  —  h)  <  L0  and  f[  and  C  are  vectors  indicating  the  location  of  the 
changes  in  Cl .  Clearly,  in  this  assumed  structure  function  of  the  refractive  index 
fluctuations,  significant  changes  in  C‘l  only  occur  over  distances  greater  than  or  equal 
to  Lq.  Assuming  that  the  atmosphere  can  be  divided  into  stacking  layers  of  thickness 
Ax:  along  the  propagation  path  such  that  the  fluctuations  in  log  amplitude  and  phase 
in  the  different  layers  are  uncorrelated,  the  wave  structure  function  after  passing 
through  N  layers  may  be  expressed  as, 

N 

D(r)  =  Y,Di(r)  (A.  100) 

i= 1 

where  D,  is  associated  with  the  ith  layer.  Taking  Zi  to  represent  the  z  coordinate  in 
the  middle  of  the  ith  layer  and  using  Equation  A. 97  for  each  individual  layer  wave 
structure  function,  the  total  wave  structure  function  may  be  expressed  as 

N 

D(r)  =  2.91k2  Y,Cl(zi)Azr*. 

i= 1 
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Generalizing  because  variations  of  C2  are  assumed  to  be  slow  compared  to  the  length 
Az,  replace  the  hnite  sum  with  the  integral  along  the  propagation  path, 

Z 

D(r)  =  2.91  k2ri  J  Cl(z')dz'  (A.102) 

o 

where  z!  is  the  variable  along  the  propagation  path  from  0  to  z.  Using  Equation  A.102 
in  the  expression  for  the  long  exposure  OTF,  Equation  A. 88,  the  long-exposure  OTF 
with  a  varying  C2  is 

Jciiz'W 

OTFL(fs)  —  exp  — 5.827T2//- - , -  .  (A.103) 

A3 

See  Section  III. A  for  details  on  the  approximations  for  how  C 2  varies  in  the 
vertical. 
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APPENDIX  B.  EXACT  FORM  EYE/BRAIN 
SPATIAL  INTEGRATION 


The  complete  expressions  for  the  eye/brain  system  spatial  integrations  are 
given  below,  assuming  the  eye/brain  spatially  integrates  over  all  frequencies  within 
the  yet  undegraded  image  of  the  target.  The  first  subscript  in  these  equations  is  the 
direction  of  spatial  integration.  The  second  is  the  orientation  of  the  MRTD  target. 

In  scanning  systems,  for  horizontal  MRTD  target  patterns,  the  system  spatial 
integrals  are  given  in  the  horizontal  and  vertical  spatial  integration  directions.  In  the 
horizontal, 


EkAf.) 


^7-  /  SMI/WOl2 

Jp  Q 


2/s 


(B.l) 


and  in  the  vertical, 


—  OO 


2/s 


(B.2) 


where  HNFhiy)  and  Hnfv{v)  are  the  horizontal  and  vertical  system  noise  filter  MTF 
components  as  given  in  Equation  2.20,  A fp  is  defined  by  Equation  2.24,  is  the 
samples  per  detector  vertical  IFOV  and  8V  is  the  detector  vertical  IFOV  in  mrad. 

For  vertical  MRTD  target  patterns  in  scanning  systems,  the  system  spatial 
integrals  are  given  in  the  horizontal  and  vertical  spatial  integration  directions.  In  the 
horizontal, 


Ehv(fs) 


Vs 

A/, 


OO 

f  S(V)[HNFh(v)f 

1 

B' 

toH 

_ _ i 

7l TV 

J 

0 

2/s 

dv 
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and  in  the  vertical, 
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For  horizontal  MRTD  target  patterns  in  staring  systems,  the  system  spatial 
integrals  are  given  in  the  horizontal  and  vertical  spatial  integration  directions.  In  the 
horizontal, 


£*„(/«) 


(B.5) 


and  in  the  vertical, 


—  OO 


2  fs 
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where  5h  is  the  detector  horizontal  IFOV  in  mrad  and  Sh  is  the  samples  per  detector 
horizontal  IFOV. 

For  vertical  MRTD  target  patterns  in  staring  systems,  the  system  spatial  in¬ 
tegrals  are  given  in  the  horizontal  and  vertical  spatial  integration  directions.  In  the 
horizontal, 


ekU.) 


2f, 


(B.7) 


and  in  the  vertical, 


—  OO 


2  fs 
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APPENDIX  C.  TEMPERATURE  SCALING  OF 

O'TVH 

In  order  to  scale  atvh,  the  model  scales  the  detector  spectral  detectivity  from 
300  K  to  T(  where  T\  is  the  new  temperature.  This  scaling  can  only  occur  if  the 
detector  noise  is  limited  by  the  incident  background  photon  flux,  in  other  words,  the 
background  limited  in  performance  (BLIP)  condition.  The  BLIP  conditions  happens 
when  the  background  photon  flux  is  much  greater  than  the  signal  flux,  causing  photon 
noise  to  be  the  dominant  noise  factor  (Driggers  1999).  The  scaling  of  the  detector 
spectral  detectivity  is  given  in  Equation  C.l, 

OO 

/QA(A,300)dA 

D*(  A,  Tx)  =  D*(  A,  300)  ^ -  (C.l) 

\  /QA(A,Ti)dA 

N  0 

where  QV(X,T )  the  spectral  photon  flux  density  where  T  =  300,  Ti.  QV(X,T )  is 
defined  in  photons/s/cm3. 

For  BLIP  detectors,  using  this  scaled  D*(X,T\),  and  for  non-BLIP  detectors 
using  the  unsealed  D*,  the  appropriately  temperature-scaled  atv  is  given  as 

/r>‘(A,300)^5(A)dA 

i&tvh^T!  =  (<7tD/i)300^ - •  (C.2) 
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APPENDIX  D.  TABLES  OF  INPUTS  TO 
FLIR92,  NVTHERMIP 


The  salient  inputs  to  NVThermlP  are  listed  in  Tables  XVI,  XVII,  and  XVIII. 
Inputs  to  NVThermlP  are  made  in  the  form  of  a  text  hie,  so  in  the  text  hie  each  entry 
is  of  the  form  ‘variable, value’  as  in  ‘Sensor  Name,  Sensor  Test  1’  if  SensorName  is  the 
variable  and  Sensor Testl  is  the  value,  or  the  name  of  the  sensor  being  tested.  Inputs 
to  FLIR92  are  listed  in  Tables  XIX  and  XX.  These  inputs  are  made  in  a  MatLab  hie. 
The  values  given  are  generic  and  correspond  to  no  special  case.  The  H  stands  for  the 
horizontal  direction  and  the  V  for  the  vertical  direction.  Atmospheric  data  changed 
for  each  case  and  is  listed  in  a  separate  Tables  X  and  XI  given  in  Chapter  IV.  Note 
that  calculations  were  carried  out  for  a  staring  sensor. 


Table  XVI.  Input  Parameters  NVThermlP  Part  I. 


Variable 

Value 

Units 

Comments 

SensorName 

Sensor  Testl 

- 

Any  name  okay 

TypeOflmager 

0 

- 

0  for  Staring 

SingleFrame 

1 

- 

- 

Cut  OnWavelcngth 

3.3 

/am 

- 

Cut  OffWavelength 

5 

yUlU 

- 

Magnihcationln 

0 

- 

0,  Calculated  in  NVThermlP 

HorizontalFOV 

0.48 

deg 

Same  as  FOV  calc,  in  FLIR92 

VerticalFOV 

0.48 

deg 

Same  as  FOV  calc,  in  FLIR92 

FrameRate 

30 

Hz 

— 

Verticallnterlace 

2 

- 

- 

Electroniclnterlace 

1 

- 

DiffWavelength 

4.15 

/jrn 

- 

AveOpticTrans 

1 

frac 

- 

FocalLengthln 

45.7 

cm 

- 

FNumberln 

0 

- 

0,  Calculated  in  NVThermlP 

Apert  ureDiameterln 

15.2 

cm 

- 

OpticsBlur 

0.01 

mrad 

- 

OpticsBlurUnits 

0 

- 

0  for  object  space 

OpticsBlurType 

0 

- 

0  for  RMS 

VibBlurX 

0.02 

mrad 

— 
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Table  XVII.  Input 


Variable 

Value 

VibBlurY 

0.02 

VibBlurType 

0 

DetHorrDim 

23 

DetVertDim 

23 

PeakDstarln 

1E+12 

IntegrationTimeln 

7000 

NumTDI 

1 

ScanEfficiency 

1 

NumHorrDet 

512 

NumVertDet 

512 

SamplesPerHIFOV 

1.4 

NoiseFactorX 

0 

NoiseFactorY 

0 

Sigmavh 

0 

Sigmav 

0 

Sigmah 

0 

FPN 

0 

PtSi 

1 

EmissionCoef 

0 

BarrierHeight 

0 

Uncooled 

1 

Dither 

0 

LowPass3dB  Cutoff 

22000 

LowPassF  ilterOrder 

1 

NoiseGain 

0 

Framelntegration 

1 

Interph 

0 

Interpv 

0 

InterpTypeh 

0 

InterpTypev 

0 

EZoom 

0 

EZoomType 

0 

Boosth 

0 

Boostv 

0 

DispType 

0 

EoMUX 

1 

HorLEDSize 

0 

Vert  LED  Size 

0 

Parameters  NVThermlP  Part  II. 


Units  Comments 


mrad 


0,  Calculated 


0,  Calculated 


0,  Calculated 


0,  Calculated 


0,  Calculated 


0,  none 


1,  no 


0,  Calculated 


0,  Calculated 


1,  no 


0,  no 


in  NVThermlP 


in  NVThermlP 


in  NVThermlP 


in  NVThermlP 


in  NVThermlP 


in  NVThermlP 


in  NVThermlP 


Number  filter  poles 


0,  no  H  interopolation 


0,  no  V  interpolation 


0,  no  H  interpolation 


0,  no  V  interpolation 


0,  no  EZoom 


0,  no  EZoom 


no  H  boost 


no  V  boost 


0,  CRT 


1,  no 
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Table  XVIII.  Input  Parameters  NVThermlP  Part  III. 


Variable 

Value 

CRTType 

0 

LED  Height 

0 

LEDWidth 

0 

DispSpotHcight 

0.02 

DispSpotWidth 

0.02 

AveDispLum 

10 

MinDispLum 

0 

SceneContTemp 

10 

DispHeight 

15.24 

DispViewDist 

70 

NumEyes 

2 

AtmTransType 

0 

AerosolModel 

0 

Mo  dclEnvir  onment 

6 

TransPerKM 

0.9 

Turbulcncelndex 

IE-14 

AtmSmoke 

1 

AtmSensorAltitudeln 

1 

AtmAlpha 

1 

AtmConcLen 

1 

NumSpectDetPtsIn 

9 

SpectDetWavelengthlnl 

3.3 

SpectDetWavelengthIn2 

3.5 

SpectDetWavelengthIn3 

3.7 

SpectDetWavelengthlnd 

4.1 

SpectDetWavelengthIn5 

4.2 

SpectDetWavelengthlnb 

4.5 

SpectDetWavelengthIn7 

4.6 

SpectDetWavclengthIn8 

4.8 

SpectDetWavelengthIn9 

5 

NormDStarlnl 

0.66 

NormDStarIn2 

0.7 

NormDStarIn3 

0.74 

NormDStarIn4 

0.82 

NormDStarIn5 

0.84 

NormDStarIn6 

0.9 

NormDStarIn7 

0.92 

NormDStarIn8 

0.96 

NormDStarIn9 

1 

Units 


Comments 


0,  RMS 


cm 


cm 


cm 


cm 


ft-Lamberts 


ft-Lamberts 


K 


cm 


cm 


km 


m2/g 


m2/g 


0,  Beer’s  Law 


Not  applicable,  Beer’s  Law 


Not  applicable,  Beer  s  Law 


not  applicable,  no  smoke 


not  applicable,  no  smoke 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


Same 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 


as  FLIR92 
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Table  XIX.  Input  Parameters  for  FLIR92  Part  I. 


Variable 

Value 

Units 

Std  Dev  Optics  Blur  H  Diarn 

0.01 

mrad 

Std  Dev  Optics  Blur  V  Diarn 

0.01 

mrad 

H  Cell  Dimension 

25 

/jrn 

V  Cell  Dimension 

25 

/jrn 

Fields  Per  Frame 

2 

- 

Detectors  Disp  H  Dir 

512 

- 

Detectors  Disp  V  Dir 

512 

- 

Sample  and  Hold  Filter  Dir 

H 

- 

H  Scene  Sample  Phase  Ang 

0 

rad 

V  Scene  Sample  Phase  Ang 

0 

rad 

H  Filter  Cutoff  Freq 

0 

cycles/mrad 

V  Filter  Cutoff  Freq 

0 

cycles/ mrad 

CCD  H  Charge  Transfer  Eff 

1 

- 

CCD  V  Charge  Transfer  Eff 

1 

- 

Active  TV  Lines  Across  Disp  Ht 

512 

- 

Display  Frame  Rate 

30 

Hz 

Std  Dev  H  Rand  Image  Motion  Amp 

0.02 

mrad 

Std  Dev  V  Rand  Image  Motion  Amp 

0.02 

mrad 

Std  Dev  H  Sinsdl  Image  Motion  Amp 

0 

mrad 

Std  Dev  V  Sinsdl  Image  Motion  Amp 

0 

mrad 

Eye  Limiting  MTF 

1 

- 

Thrshld  SNR  to  Resolve  4-bar  Tgt 

2.5 

- 

Enviro  Background  Temp 

290 

K 

Det  1/f  knee  freq 

3 

f 

Band  Avg  Optics  Transmission 

1 

- 

Detect  H  Active  Dim 

23 

/jrn 

Detect  H  Active  Dim 

23 

/jrn 

Detector  Integration  Time 

0.007 

s 

Tot  Detectors  H  Direction 

512 

- 

Tot  Detectors  V  Direction 

512 

- 

H  FIR  Digital  Filter  Type 

none 

- 

V  FIR  Digital  Filter  Type 

none 

- 

Display  Height 

0.1524 

m 

Display  Luminance 

10 

ft-lamberts 

Display  Viewing  Distance 

0.7 

m 

Peak  Det  Detectivity  300K 

1  x  1012 

cmHza  /W 

Spectral  Band  Cuton  A 

3.3 

/jrn 

Spectral  Band  Cutoff  A 

5 

/jrn 

Table  XX.  Input  Parameters  for  FLIR92  Part  II. 


Variable 

Value 

Units 

Diffraction  A 

4.15 

yUlU 

H  Focal  Length 

457 

mm 

V  Focal  Length 

457 

mm 

Optical  Aperture  H  Diarn 

152 

mm 

Optical  Aperture  V  Diarn 

152 

mm 

FPA  Readout  Frame  Rate 

30 

Hz 

Pathlcngth  (sensor  alt) 

1000 

m 
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APPENDIX  E.  ACRONYMS  USED 


Table  XXI.  Acronyms  Used  in  this  Study 


Acronym 

Name 

CCD 

Charge  Coupled  Device 

CTF 

Contrast  Threshold  Function 

FLIR92 

Forward  Looking  Infrared  92 

FOV 

Field  of  View 

FPA 

Focal  Plane  Array 

IFOV 

Instantaneous  Field  of  View 

IR 

Infrared 

MDTD 

Minimum  Detectable  Temperature  Difference 

MRTD 

Minimum  Resolveable  Temperature  Difference 

MTF 

Modulation  Transfer  Function 

NETD 

Noise  Equivalent  Temperature  Difference 

NVThermlP 

Night  Vision  Thermal  and  Image  Processing 

OTF 

Optical  Transfer  Function 

PSF 

Point  Spread  Function 

SPRITE 

Signal  Processing  in  the  Element 

TAWS 

Target  Acquisition  Weapons  Software 

TTP 

Targeting  Task  Performance 
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